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This  research  mainly  deals  with  the  quantum  efficiency  enhancement  of  QWlPs 
(Quantum  Well  Infrared  Photodetectors)  by  using  the  2-dimensional  (2-D)  planar 
metal  grating  couplers.  The  motivation  is  originated  from  the  intersubband  selection 
rule  of  quantum  well,  which  requires  the  polarization  of  incident  infrared  (IR)  radi- 
ation to  have  a component  perpendicular  to  the  quantum  well  layers.  The  unique 
features  of  the  2-D  planar  metal  grating  coupler  are  that  (a)  it  effectively  scatters  the 
normal  incident  radiation  independent  of  light  polarization  and  (b)  it  can  be  easily 
fabricated  by  using  a simple  metal  lift-off  technique.  Three  grating  structures  are 
developed  to  couple  the  normal  incident  IR  light  into  the  QWlPs.  The  first  grating 
structure  is  the  2-D  square  dot  reflection  metal  grating  with  a coupling  efficiency  of 
20%.  The  coupling  efficiency  can  be  increased  by  employing  its  complementary  geom- 
etry, the  2-D  square  aperture  mesh  metal  grating.  About  70%  of  the  normal  incident 
light  is  diffracted  into  TM  waves  available  for  QWIP  absorption.  To  solve  the  corner 
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rounding  effect  observed  in  the  square  shape  grating,  a 2-D  circular  aperture  metal 
grating  coupler  is  developed,  which  has  the  same  coupling  efficiency  of  70%  and  with 
a efficiency  of  2.7  times  higher  than  that  of  the  45°  polished  edge  QWIP  with  corre- 
sponding quantum  well  parameters.  In  addition  to  the  square  symmetry  arrangement 
required  in  the  three  grating  designs,  one  more  benefit  from  the  hexagonal  symmetry 
circular  aperture  grating  is  a wider  coupling  bandwidth.  For  each  grating  structure, 
several  universal  plots  based  on  two  normalized  parameters  \jg  and  ajg  were  pre- 
sented, where  A is  the  wavelength,  g is  the  grating  period,  and  a is  the  dimension  of 
grating  shape.  The  universal  plots  provide  a convenient  way  for  grating  design. 

Modal  expansion  and  method  of  moments  are  two  basic  approaches  used  in  the 
numerical  analysis  of  grating  couplers.  The  diffracted  angle  and  diffracted  power  are 
drawn  in  universal  plots,  the  diffracted  angle  determines  the  absorption  constant, 
combining  with  the  diffracted  power  the  spectral  quantum  efficiency  of  QWIP  can  be 
calculated. 

In  addition  to  coupling  enhancement  by  grating,  the  waveguide  geometry  QWIP 
shows  a further  improvement  in  IR  absorption.  The  design  of  QWIP  waveguide  is 
carried  out  numerically.  Together  with  the  optimum  period  of  multiquantum  wells 
and  the  blocking  layers  to  reduce  the  dark  current,  a high  performance  QWIP  is 
proposed. 
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CHAPTER  1 
INTRODUCTION 

1.1  Quantum  Well  Detectors 

There  is  a great  need  for  high  resolution,  high  sensitivity  focal  plane  array  de- 
tectors in  the  long  wavelength  infrared  spectrum  of  the  8-14  fim  atmospheric  window 
for  both  the  military  and  civilian  regimes.  The  most  widely  used  infrared  detector 
materials  are  intrinsic  Si  and  HgCdTe  (MCT).  Even  with  the  rapid  growth  of  the 
silicon-based  new  technologies  and  microelectronics  industry,  it  is  generally  recog- 
nized that  for  long  wavelength  infrared  detector  applications  III-V  compound  semi- 
conductors are  superior  to  Si  [1].  This  is  due  to  different  material  properties,  such 
as  a wide  range  of  available  bandgaps,  the  direct  bandgaps  that  facilitates  efficient 
conversion  from  light  to  electricity,  and  carriers  with  very  high  mobility.  The  dif- 
ficulties in  materials  growth,  processing,  and  other  details  of  these  technologies  of 
HgCdTe  material  have  made  it  highly  desirable  to  examine  whether  performance  can 
be  improved  by  other  material  systems  [2].  Since  the  material  and  device  technolo- 
gies in  HI-V  compound  semiconductors  are  far  more  advanced  than  HgCdTe,  it  is 
natural  to  investigate  these  HI-V  compounds  for  possible  applications  in  long  wave- 
length infrared  detection.  Recently,  IH-V  quantum  well/superlattice  structures  have 
been  used  to  construct  long  wavelength  infrared  photodetectors  [3].  The  quantum 
well  is  formed  by  using  a layer  of  a narrow  bandgap  semiconductor  (such  as  GaAs) 
sandwiched  between  two  wider  bandgap  semiconductor  layers  (such  AlxGai_a;As). 
The  motion  of  the  electron  perpendicular  to  the  layers  becomes  quantized  so  that 
localized  (in  the  direction  transverse  to  the  layers)  two-dimensional  (2-D)  subbands 
of  size-quantized  states  were  formed  inside  the  quantum  well  [4-6]. 
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These  multi-quantum  well  structures  are  usually  grown  by  molecular  beam  epi- 
taxy (MBE)  technique,  which  allows  the  exact  adjustment  of  the  most  important 
device  parameters.  By  changing  the  well  width  and  barrier  height,  detectors  can 
be  designed  with  response  peaks  from  the  mid-  to  the  far-  infrared  wavelength  re- 
gions. According  to  the  electrons  transport  approaches,  two  fundamental  structures 
of  quantum  well  infrared  photodetectors  (QWlPs)  are  illustrated  in  Fig.  1.1.  For 
the  bound-to-continuous  (BTC)  QWlPs  (Fig.  1(a)),  the  electrons  photoexcited  from 
the  bound  ground  state  by  infrared  (IR)  radiation  are  transported  through  the  ex- 
tended continuum  band  above  the  barrier  under  an  applied  electric  field  [7],  while 
for  a bound- to-miniband  (BTM)  QWIP,  as  shown  in  Fig.  1(b),  these  electrons  are 
transported  through  a global  miniband  formed  by  superlattice  barriers. 

1.2  Quantum  Well  Intersubband  Absorption 

The  absorption  in  a QWIP  is  based  on  intersubband  transition  of  n-type  doped 
quantum  well  heterostructures.  Upon  illumination  of  IR  radiation,  electrons  are  ex- 
cited from  the  ground  state  into  the  excited  state  [8].  Only  in  the  excited  state, 
electrons  moving  in  the  direction  perpendicular  to  quantum  wells  becomes  possible, 
which  produce  the  electric  signal  under  an  applied  electric  field.  Therefore,  the  basic 
operation  of  QWlPs  is  photoconductive  type,  and  the  schematic  diagram  of  a detector 
pixel  is  shown  in  Figure  1.2.  Careful  optical  design  is  needed  since  one  of  the  major 
problems  encountered  in  type-I  (n-doped)  QWlPs  is  that  the  polarization  selection 
rule  for  intersubband  transition  requires  one  component  of  the  electric  field  vector 
of  the  IR  radiation  to  be  polarized  perpendicular  to  the  quantum  well  layer  planes 
[9-11].  This  results  in  the  detectors  have  no  photo-response  for  radiation  incident  in 
the  direction  normal  to  the  quantum  well  layers.  For  focal  plane  arrays  (FPAs)  ap- 
plications, as  shown  in  Fig.  1.3,  a response  to  normal  incident  radiation  is  requisite. 
For  this  reason,  special  schemes  of  coupling  to  convert  the  transverse  electromagnetic 
wave  into  a direction  suitable  for  absorption  are  necessary.  A variety  of  techniques 
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have  been  employed,  such  as  a 45°  angle  polished  facet  in  substrate  of  the  detec- 
tor sample  as  shown  in  Fig.  1.4(a),  the  illumination  at  Brewster’s  angle  orientation 
(Fig.  1.4(b)),  a prism  coupler  [12]  in  the  detector  surface  (Fig.  1.4(c)),  a lamellar 
grating  coupler  (Fig.  1.5(a)),  and  a double  periodic  metal  grating  on  the  detector 
surface  (Figs.  1.5(b),  (c)).  The  coupling  of  the  incident  intensity  to  the  intersubband 
transitions  is  too  week  to  be  used  for  sample  oriented  at  Brewster’s  angle  or  coupled 
by  prism.  The  light  coupling  through  a 45°  polished  facet  only  allows  incident  light 
to  access  one  edge  of  the  detector,  which  is  possible  to  form  only  one-dimensional 
(1-D)  linear  arrays.  For  imaging  application  in  FPAs,  it  is  necessary  to  couple  light 
uniformly  into  the  2-D  arrays  of  such  detectors.  A 1-D  linear  grating  either  by  de- 
positing metal  lines  onto  the  detector  surface  (Figs.  1.2  and  1.5(a))  or  by  chemical 
etching  the  detector  surface  into  periodic  profile  (Fig.  1.6(a))  or  followed  by  metal 
deposition  (Fig.  1.6(b))  are  potentially  applicable  to  the  fabrication  of  2-D  detector 
arrays.  However,  these  approaches  are  polarization  dependent,  only  couple  transverse 
magnetic  (TM)  component  of  the  incident  light.  To  solve  this,  a double  periodic  cross 
grating  coupler  is  used  to  effectively  couple  the  normal  incident  IR  radiation  into  the 
doped  wells  independent  of  light  polarization. 

The  detector  array  consists  of  pixels  of  QWlPs  arranged  periodically  on  the 
detector  substrate.  To  allow  normal  incidence  illumination,  the  grating  can  be  fabri- 
cated either  on  the  substrate  side  (Fig.  1.7(a))  or  on  the  top  surface  (Fig.  1.7(b))  of 
the  detectors.  It  is  worthy  of  noting  that  the  latter  hcis  the  advantage  over  the  former 
in  that  it  eliminates  the  undesirable  pixel  cross  talk  as  shown  in  Fig.  1.7. 

1.3  Grating  Coupler 

The  above  discussions  suggest  that  diffraction  grating  is  an  effective  optical  cou- 
pling means  for  fabricating  large  area  n-type  QWIP  arrays.  This  has  been  demon- 
strated in  the  early  work  [13,14]  of  intersubband  resonances  in  electron  inversion  layer 
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of  silicon.  Goossen  et  al.  [15-17]  first  employed  the  grating  to  couple  a single  quan- 
tum well  and  gave  a brief  grating  analysis.  Then,  Hasnain  et  al.  [18]  used  a etched 
triangular  grating,  Yu  and  Li  [19,20]  used  a laminal  metal  grating,  on  QWIP  cap 
layer  to  further  study  the  grating  coupling  in  multiquantum  wells.  Hasnain  obtained 
a coupling  efficiency  about  50%  to  that  of  the  45°  edge  detectors.  Planar  1-D  metal 
stripe  gratings  were  first  studied  by  Li  et  al.  [21,22],  who  used  a grating  period  ap- 
proximately equal  to  that  of  the  peak  absorption  wavelength  in  the  detector  medium 
for  many  layer  QWIP.  This  type  of  grating  is  relatively  simple  in  fabrication.  Ap- 
plying the  metal  gratings  for  light  coupling  based  on  the  excitation  and  emission  of 
surface  plasmons  to  achieve  high  responsivity  has  been  done  by  by  Kock  et  al.  [23]. 
In  this  coupling  mechanism,  the  detector  polished  backside  was  etched  into  a slot 
profile  followed  by  a 250  A gold  film  coating.  The  surface  plasmons  are  excited  at  the 
Au-air  interface  and  leak  through  the  Au  film  to  radiate  into  the  doped  quantum  in 
several  angles.  In  addition  to  these  single  period  grating,  a wider  coupling  bandwidth 
have  been  demonstrated  by  Lee  et  al.  [11]  by  a bi-periodic  grating  arrangement. 

In  the  study  of  2-D  grating,  Andersson  and  Lundqvist  [24-27]  presented  a de- 
tailed analysis  and  measurements  on  such  a grating  coupled  QWIP.  The  calculated 
optimum  period,  cavity  shape  as  well  as  the  etching  depth  is  only  applicable  in  their 
special  designed  BTC  QWIP.  They  found  that  for  the  50  period  QWIP  with  peak  ab- 
sorption wavelength  at  9.5  fim  an  optimum  grating  consists  of  a 3 //m  period,  ~ 2//m 
square  cavity  width  and  a 1 fim  etching  depth,  which  allows  about  55%  quantum 
efficiency  for  unpolarized  incident  radiation.  Although  several  complicated  equations 
were  given,  the  grating  coupling  characteristics  of  their  grating  coupler  have  never 
been  mentioned.  Recently  Levine  [28]  employed  a random  roughened  reflecting  sur- 
face to  achieve  a significant  higher  coupling  by  utilizing  many  more  passes  of  infrared 
light  through  the  detector  active  region.  This  idea  is  originated  from  the  total  inter- 
nal reflection  of  light  trapping  on  crystalline  silicon  by  Yablonovitch  and  Cody  [29] 
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based  on  a statistical  mechanical  view  point. 

In  this  dissertation  we  mainly  deal  with  2-D  planar  metal  grating  for  QWIPs. 
The  individual  cells  are  equally  spaced  in  the  square  symmetrically  arranged  grating. 
According  to  its  operation,  the  2-D  square  symmetry  metal  grating  is  categorized 
into  two  regions.  The  function  of  the  grating  is  a qualitatively  different  one  above 
and  below  \ = g.  In  the  non- diffraction  region,  light  wavelength  ‘A’  is  greater  than 
grating  period  ^g\  the  grating  acts  as  a optical  filter  and  waves  propagate  only  in 
the  incident  and  reflected  directions,  which  are  the  zeroth  order  transmitted  waves 
and  zeroth  order  reflected  waves,  respectively  (Fig.  1.8(a)).  All  energy  that  is  not 
transmitted  by  the  grating  is  reflected  back.  In  the  diffraction  region,  light  wavelength 
‘A’  is  smaller  than  grating  period  and  the  grating  acts  as  scatters  filters.  In  this 
wavelength  region,  part  of  the  energy  may  be  diffracted  into  various  directions  of 
space  which  is  either  the  higher  order  transmitted  waves  or  higher  order  reflected 
waves  with  propagating  directions  depending  on  A as  shown  in  Fig.  1.8  (b)  and  (c). 
For  the  application  of  grating  on  QWIPs,  it  is  necessary  to  deflect  normal  incident 
IR  radiation  into  a non-zero  angle  with  respect  to  grating  normal  in  order  to  match 
the  absorption  selection  rule.  Of  course,  the  effective  components  for  absorption  are 
the  higher  order  diffracted  waves.  The  grating  diffracted  waves  consist  of  transverse 
electric  (TE)  modes  with  electric  field  vector  parallel  to  the  quantum  well  layer  planes 
as  shown  in  Fig.  1.9  (a),  and  TM  modes  with  one  component  of  the  electric  field 
vector  perpendicular  to  the  quantum  well  layers  a.s  shown  in  Fig.  1.9  (b).  Only  TM 
modes  obey  the  rule  of  intersubband  transition  and  cause  photo-signal.  Two  factors 
determine  the  strength  of  power  absorption:  (i)  cosine  of  the  angle  ‘7’  enclosed  by 
electric  field  vector  of  IR  radiation  and  the  motion  vector  of  excited  electrons,  (ii)  the 
power  density  of  effective  diffracted  waves.  It  is  preferable  to  design  the  2-D  metal 
grating  for  generating  the  higher  order  diffracted  waves  not  only  to  travel  in  a large 
cos  7 but  also  to  carry  a large  amount  of  power. 
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The  modeling  of  grating  will  be  restricted  to  thin  conducting  screen,  which  im- 
plies the  thickness  of  grating  is  much  smaller  than  light  wavelength  ‘A’.  The  first  step 
used  is  modal  expansion,  which  expands  the  unknown  electromagnetic  (EM)  waves 
near  the  grating  by  a complete  set  of  plane  waves.  It  is  also  a near-exact  approach 
for  calculating  wave  diffraction  problem,  which  considers  the  vector  characteristics 
of  the  electromagnetic  field.  Next,  method  of  moments  [30]  is  employed  to  reduce 
the  functional  equation  of  field  theory  to  matrix  equation.  Truncation  of  expanding 
series  is  made  carefully  by  examining  the  convergence  of  solution.  These  techniques 
offer  a highly  efficient  approach  for  solving  the  grating  diffraction  problem  provided 
the  grating  has  special  pixel  profiles,  e.g.,  rectangular  or  circular  shape.  With  these 
shapes,  an  analytic  expression  of  vector  function  inner  products  can  be  obtained  by 
hand  calculation. 


1.4  Synopsis  of  Chapters 

This  dissertation  contains  eight  chapters.  Chapter  1 describes  the  problem  we  are 
studying  and  reviews  ongoing  work  in  this  research.  Chapter  2 discusses  2-D  reflection 
square  dot  metal  grating.  Two  normalized  parameters  (i.e.,  s = \/g  and  h = a/g, 
where  g is  the  grating  period  and  a is  width  of  the  square  dot)  are  introduced  to 
characterize  the  2-D  reflection  metal  grating  coupler  for  QWlPs.  The  main  advantage 
for  such  a grating  structure  is  that  the  coupling  of  the  normal  incident  IR  radiation 
is  independent  of  its  polarization  direction.  Two  universal  graphs  relating  the  higher 
order  diffracted  power  and  absorption  angle  to  the  normalized  wavelength  have  been 
constructed  for  the  design  of  2-D  reflection  metal  grating  coupled  QWIPs,  which  give 
a convenient  overview  of  the  grating  characteristics  and  various  design  possibilities. 
By  using  a simple  scaling  rule,  the  characteristic  curve  of  coupling  and  subsequent 
absorption  spectrum  was  calculated. 

Chapter  3 presents  a detail  discussion  of  2-D  square  aperture  mesh  metal  grat- 
ing. Numerical  analysis  of  a 2-D  square  aperture  mesh  metal  grating  coupled  BTM 
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GaAs/AlGaAs  QWIP  has  been  carried  out.  The  main  advantages  of  this  grating 
structure  include  (1)  High  coupling  quantum  efficiency  in  the  quantum  well  absorp- 
tion spectrum,  and  (2)  coupling  of  normal  incident  IR  radiation  into  the  quantum 
wells  is  independent  of  light  polarization.  Two  normalized  parameters  s and  h (i.e., 
s = X/g  and  h = afg,  where  a is  the  width  of  the  aperture  in  the  grating)  are  used 
in  the  simulation  of  the  2-D  metal  grating  structure  formed  on  GaAs  BTM  QWIPs. 
Using  two  universal  plots  displaying  the  total  power  and  absorption  angle  of  the 
higher  order  transmission  and  reflection  diffracted  waves,  the  optical  absorption  con- 
stant and  coupling  quantum  efficiency  of  IR  radiation  can  be  calculated  for  different 
grating  periods  and  aperture  sizes  in  any  wavelength  range. 

Chapter  4 studies  the  problems  of  aperture  distortion  observed  in  the  grating 
fabrication.  The  effects  of  aperture  size  and  shape  variation  on  photoresponse  of  a 2-D 
square  mesh  metal  grating  coupler  for  a GaAs/AlGaAs  QWIP  has  been  analyzed.  The 
grating  pattern  consists  of  tiny  square  dot  distributed  on  the  detector  mesa.  During 
the  process  of  grating  pattern  photolithography,  some  distortion  in  the  aperture’s 
size  and  shape  may  occur  by  the  use  of  contact-type  mask  aligner.  This  may  lead  to 
cell-to-cell  nonuniformity  in  photo-response  for  QWIP  arrays  applications.  The  effect 
can  been  reduced  by  either  selecting  a relatively  stable  aperture  size  or  using  a larger 
grating  period. 

The  experiment  results  of  grating  coupling  were  compared  to  theoretical  calcu- 
lations in  Chapter  5.  A study  of  the  quantum  efficiency  and  responsivity  in  a GaAs 
step-BTM  QWIP  enhanced  by  the  planar  2-D  square  aperture  mesh  metal  grating  has 
be  carried  out.  Different  coupling  schemes  and  grating  parameters  were  included.  The 
responsivity  of  grating  coupled  QWIP  under  back  side  normal  incident  illumination 
is  1.65  times  higher  than  that  under  front  side.  A maximum  responsivity  of  /?/  = 0.5 
A/W  was  obtained  at  peak  response  wavelength  Ap  = 10.4//m  by  employing  a grating 
period,  g = 3.3//m  and  aperture  width,  a = 2/j,m  under  substrate  side  normal  incident 
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illumination.  The  ignoring  of  current  loss  on  the  metal  grating  surface  results  in  a 
larger  calculated  quantum  efficiency  (32%)  than  that  of  the  measured  (29%).  Com- 
bining the  measured  results  of  front  side  and  back  side  illumination  suggests  a strong 
reflection  occurring  at  the  substrate-air  interface  after  coating  the  metal  grating  on 
top  of  QWIP. 

The  coupling  of  2-D  circular  aperture  mesh  metal  grating  coupler  and  intersub- 
band absorption  of  lO^m  GaAs  BTC  QWIP  were  demonstrated  in  Chapter  6.  The 
inner  products  of  Floquet  modes  and  waveguide  modes  were  derived  in  a simple  an- 
alytic form.  Employing  normalized  wavelength  and  normalized  aperture  radius,  two 
sets  of  universal  charts  were  drawn  for  elucidating  the  normalized  total  power  and 
absorption  angle  of  the  higher-order  TM  diffracted  waves  of  the  grating.  By  the  use 
of  simple  scaling  rules,  the  optimum  grating  design  for  GaAs  QWIP  can  be  deter- 
mined. The  measured  results  show  a 2.7  times  enhancement  on  current  responsivity 
at  Ap  = lO^m  for  grating  coupled  detector  over  that  of  a 45°  angle  edge  detector.  A 
further  enhancement  can  be  expected  for  an  optimum  grating  coupled  QWIP. 

Chapter  7 develops  a standard  procedure  for  QWIP  design.  Analysis  of  the 
quantum  well  detector  parameters  for  achieving  the  optimum  performance  in  a BTM 
QWIP  was  made.  The  reduction  of  dark  current  is  achieved  by  two  simultaneous 
approaches;  a larger  barrier  height  and  a selective  blocking  layer.  To  increase  inter- 
subband absorption  in  the  quantum  wells,  a 4 fxm  AlAs  cladding  layer  beneath  the 
detector  active  region  is  employed  to  form  a waveguide  geometry  QWIP.  This  waveg- 
uide coupling  scheme  further  enhances  the  QWIP  absorption  due  to  multi-pass  of  IR 
radiation  along  the  wave  propagating  direction.  The  combination  of  1.2  ^m  cap  layer 
and  0.5  /zm  bottom  GaAs  contact  layer  ensures  the  maximum  field  intensity  occurs  in 
the  QWIP  active  region.  The  optimum  period  for  QWIPs  waveguide  deduced  from 
the  experimental  study  of  step-BTM  QWIPs  was  found  to  be  18  . Finally,  a summary 
of  this  research  is  given  in  Chapter  8. 
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Figure  1.1  Conduction  band  diagram  for  (a)  bound-to- continuous  trans 
port  QWIPs,  (b)  bound-to-miniband  transport  QWIPs. 


10 


Incident  Infrared  Light 


Figure  1.2  Schematic  representation  of  a single  detector  cell. 
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Figure  1.3  The  architecture  of  a focal  plane  arrays. 
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IR  light 


(b) 

Figure  1.4  (a)Light  coupling  of  QWIP  through  a 45°  angle  polished  facet 
on  one  edge  of  the  detector  sample,  (b)  Brewster’s  angle 
illumination  of  the  QWIPs. 
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IR  light 


Figure  1.4(c)  The  scheme  of  the  prism  coupled  quantum  well  infrared  pho- 
todetectors 
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(a) 


Figure  1.5  Three  basic  structures  of  planar  metal  grating,  (a)  1-D  strip 
metal  grating,  (b)  square  dot  reflection  metal  grating,  and 
(c)  square  aperture  mesh  metal  grating,  where  the  dark  areas 
represent  metal. 
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IR  light 


Figure  1.6  Side  view  of  the  detectors  showing  the  scheme  for  illumination 
in  normal  incidence  through  (a)  an  etched  top  contact  layer 
without  metal  deposition,  (b)  the  substrate  and  then  reflected 
by  the  etch-formed  metal  grating. 
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IR  light 


Multi-quantum  well 


(a) 


IR  light 


Figure  1.7  The  metal  grating  formed  on  (a)  substrate,  (b)  top  contact 
layer  of  the  quantum  well  detectors.  The  latter  avoids  the 
optical  crosstalk. 
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Figure  1.8  Operation  of  square  symmetry  metal  grating  in  (a)  non- 
diffraction  region,  (b)  diffraction  region  with  only  one  higher 
order  diffracted  waves,  (c)  diffraction  region  with  two  higher 
order  diffracted  waves  under  the  normal  incident  irradiation 
at  a fixed  wavelength  A. 
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Metal  grating 
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Figure  1.9  Cross  section  view  showing  the  directions  of  grating  diffracted 
waves  for  (a)  TE  diffracted  waves  with  no  electric  field  compo- 
nent in  the  carrier  motion  direction,  (b)  TM  diffracted  waves 
with  an  electric  field  component  along  the  direction  of  electron 
motion. 


CHAPTER  2 

A NUMERICAL  ANALYSIS  OF  REFLECTION  SQUARE 
DOT  METAL  GRATING  FOR  MULTIQUANTUM  WELL 
INFRARED  PHOTODETECTORS 

2.1  Introduction 

Metal  grating  coupler  operating  in  the  diffraction  region  (i.e.,  g > X ) has  drawn 
considerable  interest  in  recent  years  since  it  can  significantly  improve  the  intersubband 
absorption  in  a multiple  quantum  well  infrared  photodetector  (QWIP)  [15-18]  under 
normal  incident  illumination.  A wide  variety  of  QWIPs  relying  on  intersubband 
transitions  in  the  doped  quantum  wells  have  been  demonstrated  for  operating  in  the 
atmospheric  window  of  8 to  14  fim  spectral  range  at  77  K [31-35]. 

The  main  problem  encountered  in  a type-I  QWIP  is  that  the  electric  field  vector 
of  the  incident  infrared  (IR)  radiation  must  have  a component  perpendicular  to  the 
quantum  well  (QW)  layer  planes  to  satisfy  the  selection  rules  of  quantum  mechanics 
intersubband  absorption  [36-41].  Therefore,  the  angle  of  incident  light  with  respect 
to  the  QW  layers  must  be  non-zero  in  order  to  produce  photo-response  in  the  QWIP. 
One  way  to  couple  IR  radiation  into  the  QWIP  is  by  back  illumination  through  a 
45°  polished  facet  [42].  However,  this  scheme  only  allows  incident  light  access  to  one 
edge  of  the  detector,  providing  only  one  dimensional  (1-D)  linear  array.  For  imaging 
applications  in  focal  plane  arrays  (FPAs),  it  is  necessary  to  be  able  to  couple  IR 
light  uniformly  into  the  2-D  arrays  of  such  detectors.  An  alternative  way  of  light 
coupling  is  to  construct  a 1-D  lamellar  grating  to  scatter  light  into  a non-zero  angle 
either  by  etching  the  QWIPs  surface  to  deflect  light  [18]  or  by  coating  the  devices 
with  metal  on  one  side  [17,18,24]  to  reflect  light.  However,  these  approaches  are 
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polarization  sensitive,  and  only  couple  the  transverse  magnetic  (TM)  polarization 
components  of  the  incident  waves.  In  this  chapter,  we  present  a 2-D  double  periodic 
metal  grating  coupler  formed  on  one  face  of  the  QWIP  to  efficiently  couple  the  normal 
incident  IR  radiation  into  the  QWs.  This  planar  grating  structure  is  independent 
of  the  polarization  direction  of  the  incident  waves,  and  hence  is  suitable  for  FPA 
applications. 

Radiation  scattered  by  the  grating  consists  of  the  transverse  electric  (TE)  polar- 
ization with  the  electric  field  vector  parallel  to  the  QW  planes  and  the  TM  polariza- 
tion with  one  component  of  the  electric  field  vector  perpendicular  to  the  planes.  Only 
TM  polarization  leads  to  intersubband  absorption  in  the  QWs.  Within  the  QWIP 
absorption  region,  the  absorption  constant  depends  on  the  angle  between  the  electric 
field  vector  of  the  incident  waves  and  the  motion  vector  of  free  carriers  (i.e.,  the  di- 
rection perpendicular  to  the  QW  layer  planes),  which  is  defined  as  absorption  angle. 
The  relationship  is  the  cosine  of  the  angle  enclosed  by  these  two  vectors  [32,40].  For 
this  reason,  once  the  absorption  constant  of  45°  incident  angle  is  obtained,  the  power 
absorption  ratio  and  quantum  efficiency  of  the  grating  coupled  detector  can  be  de- 
termined by  simply  incorporating  the  power  and  angle  of  the  grating  scattered  light 
into  the  45°  angle  absorption  spectrum.  Therefore,  large  value  of  coupling  efficiency 
is  possible  by  a proper  selection  of  grating  parameters. 

The  analysis  presented  in  this  chapter  is  based  on  a planar  double  periodic  re- 
flection grating  with  metal  dots  in  a square  shape  formed  on  top  surface  of  a QWIP 
device.  Here  a normal  incident  electromagnetic  (EM)  wave  impinging  on  the  sub- 
strate side  of  the  detector  is  considered.  For  simplicity,  the  power  of  the  incident 
waves  is  assumed  to  be  equal  to  unity,  and  the  power  of  the  higher  order  diffracted 
waves  depends  only  on  the  ‘normalized  wavelength’ s (where  s = A/ff,  g is  the  grating 
period  and  A is  the  wavelength)  and  ‘strip  factor’  h (where  h = afg,  a is  the  width 
of  the  square  dot).  The  absorption  angle  used  to  determine  absorption  constant  also 
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depends  on  s.  As  a result,  to  design  such  a grating  structure  the  number  of  param- 
eters can  be  greatly  reduced.  Many  valuable  papers  have  been  published  that  are 
concerned  with  the  scattering  of  a plane  wave  by  a cross  grating  [43-46],  where  the 
reflectivity  of  the  grating  is  available,  however,  comparable  normalized  charts  relating 
the  higher  order  diffracted  mode  do  not  exit.  The  purpose  of  this  chapter  is  to  de- 
velop universal  plots  which  relate  the  power  of  higher  order  diffracted  waves  and  the 
diffracted  angle  to  the  normalized  wavelength  for  a 2-D  reflection  square  dot  metal 
grating  coupler.  Using  these  universal  plots,  optimum  grating  sizes  can  be  obtained 
for  the  2-D  reflection  metal  grating  coupler  in  any  spectral  range. 

2.2  The  Theoretical  Formalism 

The  diffracted  EM  waves  in  the  reflection  metal  grating  are  modeled  by  using 
modal  expansion  technique  [15,47].  The  grating  consists  of  arrays  of  perfect  con- 
ducting square  dots  with  infinitesimal  thickness  compared  to  the  wavelength.  These 
conducting  dots  are  arranged  periodically  along  two  orthogonal  coordinates  x and  y 
as  shown  in  Fig.  2.1.  In  this  figure  we  choose  spherical  coordinates;  6 is  the  angle 
between  the  wave  propagation  vector  k and  the  normal  to  the  plane  of  the  grating, 
and  is  the  angle  between  the  x-axis  and  the  projection  of  k on  the  x-y  plane.  The 
distribution  of  electromagnetic  field  near  the  array  is  in  the  form  of  Floquet  mode 
functions  according  to  wave  theory  in  a periodic  structure.  Under  normal  incidence, 
the  solution  of  scalar  wave  equation  is  given  by 


(2.1) 


where  the  time  dependence  is  omitted  and  the  wave  vectors  k^,  k^,  and  in 
the  X-,  y-  and  2:-  axes  are  given,  respectively,  by 
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fcf  = 


for  P > (kl)'^  + {k^y 
for  k^  < {klf  + {klf 


where  k is  the  wavenumber  in  the  incident  side,  <j>pq  and  Opq  are  the  corresponding 
spherical  coordinates  for  (p,q)  order  Floquet  mode,  and  p,q  = 0,  ±1,  ±2, . . . , ±oo. 
The  2-direction  wave  vector  k^^  depends  on  p and  q,  which  is  real  for  propagating 
modes,  pure  imaginary  for  evanescent  modes.  The  orthonormal  mode  functions  for 
the  TE  and  TM  Floquet  modes  transverse  with  respect  to  z can  be  expressed  as 
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= -{cos  (f)pqX  + sin  (f)pqy)U, 
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pq 


for  TE  modes 
for  TM  modes 
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where  g is  the  grating  period  in  both  directions  and  thus  served  as  a normalized 
factor.  The  wave  impedances  in  the  2-direction  are 
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where  Zo  = \J is  the  characteristic  impedance.  A plane  wave  with  unit  electric 
field  intensity  normal  incidence  in  the  (f>  plane  can  be  expressed  as  the  sum  of  TE  and 
TM  plane  waves,  that  is, 

2 

E'  = AoOr^OOr  (2-6) 

r=l 

where  Aoor  is  the  magnitude  of  incident  field  component  which  depends  on  the  po- 
larization direction.  The  third  subscript  r = 1 or  2 is  used  to  designate,  respectively, 
the  TE  and  TM  Floquet  modes.  Similarly,  the  scattered  field  can  also  be  expressed 
in  terms  of  the  Floquet  modes  and  reflection  coefficients  Rpqr  as  follows: 
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E = Rpqr^pqr 
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The  boundary  condition  needs  that  the  electric  field  on  the  conducting  plate  equals 
zero 

E'  + E“  = 0 at  z = 0 (2.8) 


It  has  been  shown  that  the  current  modes  nmn/  other  than  the  Floquet  modes 
were  used  to  expand  the  induced  current  —zxH’on  the  conducting  dot  in 
order  to  provide  a faster  convergence  [47].  The  limn/  are  the  dual  field  functions  of 
the  transverse  electric  field  functions  of  the  aperture  spanned  by  the  conducting  dot, 
which  satisfy  the  appropriate  boundary  conditions  on  the  dot.  For  TE  current  modes 
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where 
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is  the  normalized  factor  of  the  vector  function  Ilmni)  ^On  is  the  Newmann  factor,  Con 
equals  1 for  n = 0 and  Co„  equals  2 for  n > 1,  and  a is  the  width  of  square  conducting 
dot.  The  TM  current  modes  Yimn2  can  be  derived  in  a similar  way  as  Eq.  (2.9). 
Employing  ‘method  of  moments’,  an  equation  in  the  following  matrix  expression  can 
be  obtained 
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where  the  unknowns  Wmnl  are  the  coefficients  given  in  the  current  expression,  the 
terms 
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are  the  inner  product  of  vector  functions,  and 


Imni  — ^OOiC'toi**  + ^002  C*! 


^002 


(2.12) 


is  the  matrix  depending  on  incident  polarization.  The  reflection  coefficients  are  given 
by  substituting  current  modes  expansion  —z  x = Xlm  n X)/ into  Eq. 
(2.7).  The  result  is 

2 

Hv,r  = - W E Z E (2-13) 


m n /=1 

The  normalized  power  carried  by  (p,  q)  order  TM  propagating  mode  is  written  by 


•-’■111 


COS  ^pq2 

Since  Wmnl  varies  with  and  (7™"^  is  proportional  to  1/A^  for  a given  g,  the  reflection 
coefficient  Rpq^  in  Eq.  (2.13)  is  a function  of  \/g.  As  a result,  the  normalized  power 
Vpq2  is  also  a function  of  \/g.  The  absorption  angle  between  the  electric  held  vector 
of  the  higher  order  TM  diffracted  mode  Rpq2^pq2  and  the  z-axis  is  designated  as 
7p,2?  and  cosine  of  this  angle  depends  on  the  order  of  diffraction  and  the  normalized 
wavelength  s{—  X/g^  where  A is  the  wavelength  in  the  detector  medium)  as  given  by 


cos  7p,2  = S^p2  -t- 


(2.15) 


2.3  Results  and  Discussion 


As  discussed  above,  the  evanescent  modes  of  the  grating  do  not  produce  intersub- 
band absorption  [17]  in  the  quantum  wells;  neither  do  the  TE  modes.  Therefore,  the 
following  discussion  will  focus  only  on  the  propagating  TM  modes.  It  is  convenient 
to  introduce  two  normalized  parameters  here  for  use  in  the  universal  plots  for  the 
square  dot  metal  grating  structure;  these  are  the  normalized  wavelength  s = Xjg  and 
the  strip  factor  h = a/g.  In  the  computer  implementations,  for  reasons  of  available 
memory  space  (limited  to  about  20,000  double  precision  complex  numbers),  the  series 
in  Eqs.  (2.13)  and  (2.7)  were  truncated  so  that  up  to  forty  current  modes  and  5000 


25 


Floquet  inodes  were  included.  The  numerical  results  presented  below  are  sufficiently 
accurate  by  checking  their  convergence  with  increasing  numbers  of  Floquet  modes 
and  current  modes.  For  larger  strip  factors  {ajg  > 0.85  ),  a much  larger  number  of 
current  modes  is  required  to  expand  the  induced  current.  Analysis  of  large  strip  fac- 
tor are  also  unnecessary  since  it  is  more  difficult  to  fabricate  such  grating  on  QWIP. 
The  incident  waves  are  plane  waves  that  have  an  angle  of  zero  degree  with  respect 
to  the  array  normal.  In  order  to  interact  with  the  doped  quantum  well,  the  effective 
coupling  of  the  diffracted  wave  is  due  to  TM  Floquet  modes.  For  grating  period  g 
less  than  wavelength  A,  the  grating  is  operated  in  non-diffraction  region  with  only 

—4 

one  (zeroth  order)  far  field  propagating  mode  Rooi^ooi,  i = 1,2.  The  approach  used 
here  to  create  effective  coupling  is  to  excite  the  higher  order  TM  diffracted  waves  into 
the  QWs.  At  A = 5^  the  four  (1,0)  diffracted  waves  Rou^o^,  ^-12^0-^,  R-i02^^w 
and  i?io2^io^  become  propagating.  At  longer  wavelength  they  are  evanescent  and 
their  amplitude  is  attenuated  to  1/e  at  a distance  {g — g^  j X^)~  on  the  grating 
dot  matrix.  For  g ^ X,  this  distance  is  larger  than  the  spacings  of  the  dot,  and 
the  evanescent  waves  can  transfer  substantial  amounts  of  energy  from  one  dot  to  the 
next.  Therefore,  a,t  X = g the  total  power  carried  by  propagating  waves  drops  to  a 
minimum  value.  Figure  2.2(a)  shows  the  first  universal  plot,  which  illustrates  that 
the  total  power  of  the  first  order  TM  diffracted  waves  as  a function  of  the  normalized 
wavelength  for  s < 1.  For  a square  conducting  dot  grating,  the  total  power  of  the  first 
order  TM  diffracted  waves  varies  with  s = Xfg  but  is  independent  of  the  polarization 
direction  of  the  incident  waves.  This  is  due  to  the  structure  symmetry  in  the  x and 
y directions  under  normal  incidence.  Figure  2.2(b)  shows  a plot  of  total  diffracted 
power  versus  wavelength  for  a grating  with  g = 4/im  for  different  dot  width  a. 

The  absorption  constant  depends  essentially  on  the  angle  between  the  electric 
field  vector  and  the  normal  of  the  QW  layers  within  the  absorption  band  [32,40]. 
For  a given  QW  structure,  the  total  power  absorbed  by  the  detector  can  be  easily 
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determined.  Figure  2.2(a)  shows  that  the  total  power  of  the  first  order  TM  diffracted 
waves  will  increase  with  increasing  the  value  of  h.  Figure  2.3  illustrates  the  total 
power  of  the  second  order  TM  diffracted  components  and  i?i±i2^i±i2 

as  a function  of  the  normalized  wavelength  s for  different  values  of  h.  A comparison  of 
Figs.  2.2(a)  and  2.3  reveals  that  the  first  order  diffraction  power  is  much  larger  than 
the  second  order  diffraction  power,  particularly  for  larger  values  of  h.  The  second 
order  diffracted  modes  emerge  for  s < 1 / \/2  with  total  power  about  one  eighth  that 
of  the  first  order  diffracted  modes. 

Figure  2.4  shows  an  universal  plot  that  relates  cos7p,2  to  the  normalized  wave- 
length of  the  higher  order  diffracted  waves,  where  7 is  the  angle  enclosed  by  the 
electric  field  vector  of  the  diffracted  waves  and  the  grating  normal.  The  relationship 
for  the  first  order  diffracted  waves  ( i.e.  | p | -f  | 9 |=  1 in  the  notation  of  Floquet 
modes  ) is  given  by 


cos7p,2  = - (2.16) 

9 

The  relationship  for  the  second  order  diffracted  waves  with  | p | = | 9 |=  1 is 

cos7p,2  = V^s  (2-17) 

where  A is  the  wavelength  in  the  detector  medium. 

Due  to  the  dependence  of  the  absorption  constant  on  cos^7p,2,  we  suggest  that 
the  metal  grating  period  to  be  selected  in  the  region  where  s varying  between  0.7 
and  1 when  utilizes  the  first  order  diffracted  waves.  Referring  to  Fig.  2.2(a),  it  is 
noted  that  for  /i  = 0.8  the  coupling  is  relatively  flat  over  a broad  wavelength  regime. 
Thus,  h = 0.8  may  be  considered  as  a better  selection  for  coupling  IR  light  into  the 
QWs.  For  example,  for  a GaAs  QWIP  with  peak  response  Ap  at  10  pm  wavelength, 
by  choosing  s = Ap/p  = 0.77  and  h = 0.8,  we  obtain  a relatively  flat  effective  coupling 
spectrum  covering  the  wavelength  from  8 to  12  pm  as  shown  in  Fig.  2.2(b).  Since  the 
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refraction  index  of  GaAs  contact  layer  is  equal  to  3.25  at  77/t',  the  grating  period 
and  metal  dot  width  are  found  to  be  = 4.0/zm  and  a = 3.2pm,  respectively. 

According  to  Refs.  49  and  50,  the  absorption  constant  can  be  fit  by  a Lorentzian 
line  shape  in  the  following  expression 


or(i/,7) 


oc 


cos^  7 


(2.18) 


[{u-uoY  + T^t^)] 

where  r(i/)  is  the  line  width,  uq  is  the  peak  absorption  frequency  and  7 is  the  ab- 
sorption angle.  The  T(i/)  and  i/q  were  used  as  two  fitting  parameters.  The  QWIP 
structure  used  here  as  an  example  of  grating  design  is  the  extended  long  wavelength 
A = 11  ~ 15pm  GaAs/AlGaAs  QWIP  based  on  a localized  bound  state  -to-  extended 
continuous  band  transition  reported  by  Zussman  et  ah  [50],  which  consists  of  50 
periods  of  GaAs  quantum  wells  with  well  width  L-u,  = 60A  (doped  n = 5 x 10^^ 
cm~^)  and  Alo.15Gao.85As  barriers  of  thickness  Lb  = 500A  (i.e.,  a superlattice  length 
I = 2.75pm),  sandwiched  between  0.5pm  top  and  1pm  bottom  contact  layers  GaAs 
of  n = 1 X 10^*  cm~^  GaAs.  After  comparing  the  45°  incident  absorption  constant 
given  by  Ref.  51,  the  Lorentzian  fit  is  plotted  in  Fig.  2.5  as  a broken  line.  In  the 
figure,  we  also  plot  the  absorption  constant  of  the  above  detector  sample  coupled  by 
a square  dot  metal  grating  under  normal  incident  illumination,  where  g = 5pm  is  the 
grating  period  and  a = 4pm  is  the  square  metal  dot  width.  The  quantum  efficiency 
[51]  corresponding  to  the  absorption  curve  of  Fig.  2.5  is  illustrated  in  Fig.  2.6.  The 
low  quantum  efficiency  found  in  Fig.  2.6  is  mainly  due  to  the  low  doping  density  of 
the  QWIP. 


The  recent  results  on  a grating  coupled  QWIP  shows  a good  agreement  with 
the  theoretical  analysis.  The  QWIP  used  in  our  experiment  was  an  InGaAs/InAlAs 
miniband  transport  QWIP  [52]  coupled  by  a square  dot  reflection  metal  grating  with 
grating  parameters  g = 4pm  and  a = 3pm  under  normal  incidence  illumination.  The 
refractive  index  of  the  detector  medium  is  equal  to  3.7.  A comparison  between  the 
theoretical  calculations  and  the  experimental  data  is  illustrated  in  Fig.  2.7. 
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2.4  Conclusion 

In  conclusion,  we  have  performed  a numerical  analysis  of  the  normal  incident 
2-D  reflection  square  dot  metal  grating  coupler  for  the  QWIP  devices.  The  method 
of  moments  is  used  to  implement  the  numerical  calculations.  The  main  advantage 
for  such  a grating  structure  is  that  the  coupling  of  the  normal  incident  IR  radiation 
is  independent  of  light  polarization  direction.  It  is  shown  that  the  total  power  and 
the  angle  of  the  higher  order  diffracted  waves  depend  on  two  normalized  parameters 
s = Xfg  and  h = afg  which  are  functions  of  wavelength  and  square  dot  dimensions. 
By  using  simple  scaling  rule  in  the  universal  plots  shown  in  Figs.  2.2  to  2.4,  the 
optimized  grating  period  and  metal  dot  width  can  be  obtained  for  a specified  infrared 
band  detection  in  imaging  applications. 
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THIN  CONDUCTING  DOTS 


(a) 


Figure  2.1  Schematic  diagram  showing  the  square  dot  reflection  metal 
grating  and  the  directions  of  incident  and  refracted  waves, 
(a)  Top  view,  (b)  Side  view. 
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2.22.0  1.8  1.6  1.4  1.2  1.0 


normalized  wavelength  s = ^/g 


Figure  2.2(a)  Total  power  of  the  first  order  diffracted  waves  versus  s = X/g 
for  different  values  of  h = a/ g,  where  g is  grating  period  and 
a is  metal  dot  width. 
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free  space  wavelength  X (Mm) 


Figure  2.2(b)  Replot  of  Fig.  2.2(a)  versus  wavelength  for  g = 4.0  gm  and 
a = 2.0, 3.2  gm.  The  refractive  index  of  GaAs  is  3.25  in 
IIK.  Normal  incident  backside  illumination  of  the  grating  is 
assumed. 
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g/x 

2.2  2.0  1.8  1.6  1.4  1.2  1.0 


normalized  wavelength  s = ^/g 


Figure  2.3  Total  power  of  the  second  order  diffracted  wave  as  a function 
of  s = \jg  for  different  values  of  h = a/g.  Normal  incident 
backside  illumination  of  the  grating  is  assumed. 
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normalized  wavelength  s = ^/g 


Figure  2.4  Cosine  of  the  angle  between  the  electric  field  vector  of  TM 
diffracted  waves  and  the  z-axis  versus  normalized  wavelength. 
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wavelength  X {jjm) 


Figure  2.5  The  absorption  constant  of  an  extended  long  wavelength 
QWIP  [50]  by  both  grating  coupling  and  45°  light  incidence. 
Where  g = 5/im  is  the  grating  period  and  a = 4/im  is  the 
width  of  square  dot  in  the  metal  grating. 
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Figure  2.6  The  corresponding  coupling  quantum  efficiency  of  Fig.  2.5. 
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Figure  2.7  The  responsivity  of  a 2-D  square  dot  metal  grating  coupled 
BTM  QWIP  measured  at  61  K.  The  QWIP  consists  of  20 
periods  of  enlarged  Ino.53Gao.47As  quantum  wells  with  a well 
width  of  IlOA  and  a dopant  density  5x10^^  cm“^.  The  bar- 
rier layers  on  each  side  of  the  QW  consist  of  6 periods  of 
undoped  Ino.52Alo.48As(35A)/Ino.53Gao.47As(50A)  superlattice 
layers.  The  grating  period  g = 4//m  and  dot  width  a = 3/xm. 


CHAPTER  3 

DESIGN  OF  A TWO-DIMENSIONAL  SQUARE  APERTURE  MESH 
METAL  GRATING  COUPLER  FOR  A MINIBAND  TRANSPORT 
GaAs  QUANTUM  WELL  INFRARED  PHOTODETECTOR 

3.1  Introduction 

Recently,  a considerable  interest  in  the  long  wavelength  intersubband  quantum 
well  infrared  photodetectors  (QWIPs)  for  operating  in  the  8 to  12/im  atmospheric 
window  region  at  77K  has  been  reported  in  the  literature  [53-56].  Most  of  these 
QWIPs  were  fabricated  using  a large  bandgap  HI-V  system  semiconductor  material 
such  as  GaAs/AlGaAs  instead  of  the  more  difficult  narrow  bandgap  material  such  as 
HgCdTe.  The  GaAs  QWIP  offers  a very  promising  approach  for  the  long  wavelength 
infrared  (IR)  detection  due  to  the  matured  GaAs  growth  and  processing  technology. 
Low  cost  and  extremely  uniform  large  area  focal  plane  arrays  (FPAs)  can  be  fabricated 
using  GaAs  QWIPs  for  staring  IR  image  sensor  applications  [57].  However,  these 
intersubband  transition  QWIPs  do  not  absorb  normal  incident  IR  radiation,  since 
the  electric  field  vector  of  the  incident  light  must  have  a component  perpendicular  to 
the  quantum  well  (QW)  layers  in  order  to  induce  intersubband  transitions  [39,58-59]. 
As  a result,  the  angle  of  incidence  with  respect  to  the  normal  of  QW  layers  must 
be  different  from  zero  in  order  to  produce  photo-signal.  On  the  other  hand,  for  the 
fabrication  of  FPAs,  a response  in  normal  incident  light  is  required.  A linear  grating 
is  candidate  for  applicant  to  2-D  detector  array  [17,21,22],  nevertheless,  it  suffers 
from  a low  coupling  efficiency  due  to  its  polarization  selectivity.  A double  periodic 
metal  grating  formed  on  the  top  of  a QWIP  can  be  employed  to  deflect  the  normal 
incident  light  into  an  absorbable  angle  independent  of  light  polarization.  There  are 
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two  approaches  to  form  a double  periodic  grating  on  a QWIP:  one  method  is  to  deposit 
a metal  grating  directly  onto  the  detector  as  shown  in  Fig.  3.1,  and  the  other  approach 
is  to  chemically  etch  the  top  surface  layer  of  the  detectors  to  a desired  profile  followed 
by  metal  deposition  [25,26].  The  former  has  the  advantage  over  the  latter  in  that  it 
uses  a simple  planar  metal  grating  structure  which  avoids  the  unwanted  nonuniform 
etching  on  the  detector  surface  and  hence  reduces  the  processing  difficulties.  In  this 
chapter  we  report  a numerical  analysis  of  a planar  2-D  square  mesh  metal  grating  by 
depositing  metal  grid  directly  onto  the  top  cap  layer  of  the  GaAs  bound-to-miniband 
(BTM)  QWIP.  The  detector  can  be  illuminated  either  on  the  front  side  or  the  back 
side.  The  results  showed  that  coupling  quantum  efficiency  of  the  detector  under  back 
side  illumination  is  2.34  times  higher  than  that  under  the  front  illumination. 

For  practical  applications,  the  electromagnetic  (EM)  waves  impinge  on  a QWIP 
under  normal  incident  illumination.  The  radiation  scattered  by  the  2-D  metal  grating 
consists  of  the  TE  polarization  with  electric  field  vector  parallel  to  the  QW  planes  and 
the  TM  polarization  with  one  component  of  the  electric  field  vector  perpendicular  to 
the  QW  planes.  Only  TM  components  of  the  IR  radiation  will  lead  to  intersubband 
absorption  in  the  QWs.  In  the  square  mesh  metal  grating  structure  shown  in  Fig.  3.1, 
for  a unit  amplitude  normal  incident  light,  the  total  power  of  each  order  diffracted 
mode  depends  on  the  ‘normalized  wavelength’ s = \ fg  (where  g is  the  grating  period 
and  A is  the  free  space  wavelength)  and  the  ‘strip  factor’  h = a/ g (where  a is  the 
width  of  the  square  aperture  in  the  metal  grating).  Furthermore,  the  angle  between 
each  order  of  the  diffracted  waves  and  the  grating  normal  defined  by  absorption 
angle  is  determined  by  the  parameter  s and  the  diffracted  order.  The  intersubband 
absorption  in  the  quantum  well  of  a BTM  QWIP  is  a function  of  the  total  power  of  the 
higher  order  TM  diffracted  waves  and  the  diffracted  angle  [32].  The  grating  scattering 
problem  has  already  been  investigated  by  many  authors  [43-46].  The  objective  of  this 
chapter  is,  however,  not  to  present  yet  another  solution  to  this  scattering  problem 
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but  to  generate  the  novel  universal  plots  from  which  the  total  power  of  higher  order 
TM  diffracted  waves  and  the  absorption  angle  can  be  determined  for  different  grating 
periods  and  aperture  sizes  of  the  2-D  square  aperture  mesh  metal  grating  coupler  on 
GaAs  BTM  QWIP.  Together  with  the  absorption  constant  of  intersubband  transition, 
these  plots  provide  a convenient  way  for  obtaining  an  optimum  grating  coupler  for 
the  BTM  QWIPs  structure. 


3.2  Basic  Theory 

The  basic  theory  used  in  deriving  the  unknown  scattered  waves  of  a 2-D  square 
aperture  mesh  metal  grating  coupler  is  based  on  the  modal  expansion  technique  [59]. 
The  grating  consists  of  an  infinitesimal  thickness  perfect  conducting  screen  perforated 
with  square  apertures  distributed  periodically  along  two  orthogonal  coordinates  x and 
y as  shown  in  Fig.  3.1.  The  direction  of  scattered  waves  is  specified  by  the  spherical 
polar  angles;  6 is  the  angle  of  declination  of  the  wave  propagation  vector  k measured 
from  the  +z  axis,  and  (j)  is  the  azimuthal  angle  measured  with  a counterclockwise 
sense  from  x axis  of  Fig.  3.1(a).  The  mesh  metal  grating  is  formed  on  top  of  the 
QWIP.  The  basic  form  of  electromagnetic  field  near  the  grating  array  is  Floquet  mode 
functions,  therefore,  under  normal  incident  illumination  the  solution  of  scalar  wave 
equation  with  time  dependence  e-'"*  omitted  is  given  by 


Up,  = (3.1) 

where  the  wave  vectors  k^,  and  k^'^  in  the  x-,  y-  and  2-  axes  are  given,  respec- 
tively, by  Eqs.  (2.2)  and  (2.3)  where  k is  the  wavenumber  in  free  space.  The  vector 
orthonormal  mode  functions  for  the  TE  and  TM  modes  transverse  with  respect  to  z 
can  be  expressed  as  Eq.  (2.4).  The  wave  admittances  looking  into  the  air  region  from 
z = 0'^  plane  are 

iTE  _ ^ J_ 
k Zo 
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where  Zq  = yj Ho/^o  is  the  free  space  characteristic  impedance.  Assuming  that  the 
impedance  between  all  the  quantum  well  layers  are  perfectly  matched  and  the  inter- 
subband absorption  is  complete  so  that  no  waves  are  reflected  from  the  transmitted 
end  of  the  QWIP.  Under  this  condition,  the  modal  admittances  for  TE  and  TM  waves 
looking  into  the  QW  region  from  2 = O'*"  plane  are  obtained  by  replacing  Zq  in  Eq. 
(3.2)  by  Zd  = for  GaAs,  which  yields 

dTE  ^ ^ J_ 

k Zd 

C = 

A plane  wave  with  unit  electric  field  intensity  normal  incident  in  the  (j)  plane  can 
be  expressed  as  the  sum  of  TE  and  TM  plane  waves  as  described  in  Eq.  (2.6), 
E'  = X/r=i  ^oor^oor)  wherc  the  third  subscript  r = 1 or  2 is  used  to  designate  the  TE 
and  TM  Floquet  modes,  respectively.  In  fact,  the  wave  expansion  consists  of  both 
upward  and  downward  waves;  in  a similar  way,  the  reflected  waves  and  transmitted 
waves  can  also  be  expressed  in  terms  of  the  Floquet  modes  with  reflection  coefficient 
Rpqr  and  transmission  coefficient  Tpqr  as  following 

2 2 

E*  = ^OOr^OOr  + Rpqr^pqr 

T=l  p q r=l 

2 

= Tpqr^pqr  (3-4) 

p q r=l 

Equation  (3.4)  includes  the  boundary  condition  that  the  tangential  component  of 
electric  field  vector  in  the  aperture  is  continuous.  The  orthonormal  waveguide  modes 
Umn/  of  the  square  aperture  itself  other  than  the  Floquet  modes  were  used  to 
expand  the  unknown  electric  field  distribution  in  the  aperture  in  order  to  satisfy  the 
boundary  condition,  which  has  been  shown  to  provide  a faster  convergence  [59].  For 
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is  the  normalized  factor  of  the  vector  function  Ilmni)  and  con  is  the  Newmann  factor, 
£on  equals  1 for  n = 0 and  eon  equals  2 for  n > 1,  a and  h are  the  width  of  the  square 
aperture  in  the  direction  x and  y,  respectively,  here  b = a.  The  TM  waveguide  modes 

Hmn2  are  derived  in  a similar  way.  As  a result,  can  be  expressed  by 

2 

m n /=1 

where  the  unknowns  Wmnl  are  the  coefficients  given  in  the  waveguide  mode  expres- 
sion. After  substituting  Eq.  (3.6)  into  (3.4)  followed  by  taking  inner  product  of  vector 
functions  !!„„/  and  a matrix  equation  is  obtained  by 


L TTXYli  J 


[lEmn/]  — [Ann/] 


(3.7) 


where 
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and  the  incidence  depending  on  matrix 
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(3.9) 


The  higher  order  transmission  coefficient  and  reflection  coefficient  are  given,  respec- 
tively, by 

2 

= E E E '^"•"<'^,7'  P.  9 / 0, 0 (3.10) 

m n /=1 
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and 


2 

= p,9?^0.0  (3.11) 

m n /z=l 


Since  Wmni  is  proportional  to  and  varies  as  1/A^  for  a given  g,  the  trans- 
mission coefficient  T^qr  and  reflection  coefficient  R^qr  are  a function  of  \/g.  The  nor- 
malized diffracted  power  of  (p,  q)  order  TM  propagating  mode  is  given  as  Eq.  (2.14). 
The  absorption  angle  between  the  electric  held  of  the  higher  order  TM  diffracted 
modes  Tpq2^pq2  or  Rpq2^pq2  and  the  z direction  is  designated  as  'ypq2,  and  the  cosine 
of  this  angle  depends  on  the  diffracted  order  and  normalized  wavelength  s (=  X/g) 
as  given  by 


cos7p,2  = — \/p^  + 9^  (3.12) 

Tt‘P 

where  is  the  refractive  index  of  the  medium  in  which  the  EM  waves  propagate. 

The  2-D  square  mesh  grating  coupled  GaAs/AlGaAs  BTM  QWIP  structure  used 
in  the  present  analysis  is  based  on  the  transition  from  the  ground  bound  state  in  the 
enlarged  well  to  the  global  miniband  in  the  superlattice  barrier  [19]  as  shown  in  Fig. 
3.2.  In  this  QWIP  the  2-D  electrons  in  the  miniband  can  move  freely  in  the  direction 
perpendicular  to  the  QW  layers  and  may  give  rise  to  an  electrical  current  through 
the  QWIP.  The  intersubband  absorption  constant  can  be  expressed  by  [32] 

/ 


a = 


Tie  COS^  7 T^\J{E  - Emin){Emax  ~ E) 


m*^nrtQc)  hu;  {U  + SiY  — Si{E„ 


E) 


(3.13) 


where  e is  the  electronic  charge,  m*  is  the  effective  mass,  c is  the  speed  of  light,  Ug  is 
the  electron  density  in  the  wells,  u>  is  the  angular  frequency  of  the  EM  waves,  7 is  the 
angle  between  the  electric  held  vector  of  the  IR  radiation  and  the  motion  vector  of 
electrons,  Emax  and  Emin  are  the  two  extreme  edge  of  the  miniband  as  shown  in  Fig. 
3.2.  All  other  parameters,  Tt,U  and  Si  in  Eq.  (3.13)  depend  on  the  specific  quantum 
well  structure  and  wavelength  of  the  incident  IR  radiation.  Substituting  Eq.  (3.12) 
into  (3.13),  the  absorption  constant  of  the  QWIP  in  the  spectral  range  of  interest  can 
be  readily  calculated. 
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3.3  Results  and  Discussion 

The  evanescent  modes  excited  by  the  grating  produce  no  photoresponse  in  quan- 
tum well  intersubband  transition;  neither  do  TE  modes  [17,60].  Therefore,  the  fol- 
lowing discussion  will  be  focused  on  the  propagating  TM  modes  in  the  QW  region. 
It  is  convenient  to  use  two  normalized  parameters,  namely,  ‘s  = A/^r’  and  ‘h  = ajg' 
to  illustrate  the  universal  plots  for  the  square  aperture  mesh  metal  grating.  Rigorous 
calculations  were  made  by  using  the  lowest  40  waveguide  modes,  the  addition  of  more 
higher  order  modes  made  no  noticeable  change  in  the  transmission  and  reflection  co- 
efficients. When  analyzing  the  mesh  grating  with  larger  strip  factor  {ajg  > 0.85),  a 
much  larger  number  of  waveguide  modes  is  required  to  simulate  the  wave  distribu- 
tion at  the  metal  edge.  Large  strip  factor  is  also  undesirable  for  the  practical  QWIP, 
since  it  is  more  difficult  to  deposit  such  a fine  metal  line  (i.e.  less  than  1 //m)  on  the 
detector  surface. 

An  EM  wave  with  unit  power  density  impinging  on  the  mesh  metal  grating  at 
zero  degree  angle  with  respect  to  the  grating  normal  is  considered.  Besides,  the 
following  discussion  is  mainly  based  on  the  back  side  illumination  as  shown  in  Fig. 
3.1(b).  The  effective  coupling  is  due  to  non-zeroth  order  reflected  TM  Floquet  modes, 
since  the  TE  Floquet  modes  have  an  angle  7 = 90°  in  Eq.  (3.13).  Figure  3 shows 
the  first  universal  plot,  which  illustrates  the  normalized  total  power  of  the  first  order 
TM  diffracted  waves  Roi2^oi2>  Ro-i2^o-i2)  R— 102^— 102  3-nd  Rio2^io2  ^.s  a function 
of  the  normalized  wavelength  s = X/g  for  different  values  of  h.  We  can  also  obtain 
the  power  of  the  transmitted  Floquet  modes  under  front  illumination  shown  in  Fig. 
3.1(c)  by  dividing  the  normalized  power  in  Fig.  3.3  by  a factor  2.34,  which  includes  a 
modification  factor  0.72  to  account  for  the  power  reflection  owing  to  the  impedance 
mismatch  between  air  (n^  = 1)  and  GaAs  (nr— 3.25  at  77 K).  The  first  order  diffracted 
waves  become  propagating  when  the  wavelength  of  the  light  in  GaAs  is  smaller  than 
the  grating  period,  that  is,  s = X/g  < 3.25.  For  the  same  reason,  within  the  spectral 
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range  in  Fig.  3.3,  only  zeroth  order  far  field  transmitted  waves  Tooi^ooi  (*  = 1, 2)  in  the 
free  space  might  be  found,  and  all  other  higher  order  transmitted  waves  were  cut-off. 
This  is  due  to  the  fact  that  free  space  wavelength  A is  greater  than  the  grating  period 
g and  thus  makes  the  grating  operating  in  the  non-diffraction  region  on  the  air  side. 
In  addition,  the  square  mesh  grating  is  indistinguishable  between  x and  y directions. 
The  total  normalized  power  of  the  first  order  diffracted  waves  generated  by  x and  y 
components  of  the  incident  waves  remain  the  same  for  different  input  polarizations. 
In  other  words,  the  coupling  of  the  grating  is  polarization  independent.  A rather  flat 
coupling  curve  for  h = ajg  — 0.5  in  Fig.  3.3  may  be  treated  as  a better  choice  for 
coupling  the  IR  radiation  more  efficiently  into  the  QWIP.  Figure  3.4  illustrates  the 
normalized  total  power  of  the  second  order  TM  diffracted  components  i?±n2^±ii2 
and  -Ri±i2^i±i2  as  a function  of  the  normalized  wavelength  s for  various  values  of 
h.  Similarly,  the  second  order  transmitted  TM  components  under  front  illumination 
(Fig.  3.1(c))  is  obtained  by  dividing  the  coupling  power  in  Fig.  3.4  by  a factor  of  2.34 
if  the  reflective  effect  at  the  air-GaAs  interface  is  modified.  A comparison  of  Figs.  3.3 
and  3.4  reveals  that  those  gratings  exciting  larger  power  of  the  first  order  diffracted 
waves  will  excite  larger  power  of  the  second  order  diffracted  waves  too.  The  second 
order  diffracted  waves  emerge  for  s < 2.298  with  a total  power  about  50%  smaller 
than  that  of  the  first  order  diffracted  waves. 

Figure  3.5  shows  a universal  plot  that  relates  cos7p,2  to  the  normalized  wave- 
length of  the  higher  order  diffracted  waves,  where  7 is  the  absorption  angle.  The 

relationship  for  the  first  order  diffracted  waves  ( i.e.  | p | -(-  | 9 |=  1 in  the  notation 
— # 

of  Floquet  modes  ) is  given  by 


X Ur 

cos  7p,2  = 

9 

The  relationship  for  the  second  order  diffracted  waves  with  | p | = | ^ |=  1 is 


(3.14) 
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(3.15) 


We  next  calculate  the  absorption  constant  versus  wavelength  in  a 2-D  square 
mesh  metal  grating  coupled  GaAs  BTM  QWIP.  The  mesh  grating  period  is  selected 
in  the  region  where  2.298  < s < 3.25  falls  in  the  QWIP  absorption  band  by  using 
the  first  order  TM  diffracted  waves.  In  this  case,  cos^7  is  greater  than  1/2  which 
corresponds  to  the  light  launching  at  a 45°  angle.  If  we  select  the  curve  with  h = 0.5 
in  Fig.  3.3,  then  by  multiplying  the  grating  period  g = 4/zm  to  the  coordinate,  a 
rather  fiat  curve  of  the  coupling  efficiency  over  a broad  wavelength  regime  9-13/xm 
can  be  obtained,  as  shown  in  Fig.  3.6(a).  By  substituting  the  corresponding  cos  7 
in  Fig.  3.6(a)  and  using  the  QWIP  parameters  given  in  Fig.  3.2  into  Eq.  (3.13), 
we  obtain  the  absorption  constant  versus  wavelength  for  the  GaAs  BTM  QWIP  as 
shown  in  Fig.  3.7(a).  Note  that,  the  dot  line  represents  the  absorption  constant 
of  the  QWIP  with  IR  launched  at  a 45°  polished  facet  and  the  electron  interactive 
effect  [61,62]  has  been  considered.  The  aperture  width  of  the  mesh  metal  grating  in 
the  above  example  was  found  to  be  a =(4  /im)x0.5  = 2/xm.  The  coupling  quantum 
efficiency  g for  the  QWIP  [32,51]  can  be  calculated  by  using  the  expression 


(3.16) 


where  Veff  is  the  effective  coupling  power  of  the  light  and  / is  the  total  length  of 
the  doped  quantum  wells.  In  the  present  case,  Veff  is  the  normalized  total  power  of 
the  first  order  diffracted  waves,  I = (8SA)  x (40  periods)  = 3520A,  and  g is  shown  as 
a broken  line  in  Fig.  3.7(b). 

The  optimum  square  aperture  mesh  metal  grating  used  for  the  above  BTM  QWIP 
is  g = 3.3//m  and  a = 2/xm,  corresponding  to  the  curve  with  hjg  = 0.6  in  Fig.  3.3. 
Following  the  same  procedure,  the  characteristic  curve  of  the  first  order  TM  diffracted 
waves  for  such  a grating  under  back  illumination  was  plotted  in  Fig.  3.6(b),  which 
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shows  that  the  cut-oif  wavelength  of  the  grating  is  10.725/xm.  The  solid  line  in  Fig. 
3.7(a)  illustrates  the  absorption  constant  of  the  optimum  grating  design.  Figure  3.7(b) 
shows  the  quantum  efficiency  of  the  GaAs  BTM  QWIP  under  different  light  coupling 
schemes.  The  results  show  that  the  coupling  quantum  efficiency  of  a 2-D  square  mesh 
metal  grating  under  back  illumination  is  superior  to  that  under  front  illumination. 
Furthermore,  the  optimum  grating  has  a coupling  efficiency  larger  than  that  of  a 45° 
angle  light  coupling. 


3.4  Conclusion 

In  conclusion,  we  have  performed  a detailed  analysis  of  a GaAs  BTM  QWIP 
coupled  with  a 2-D  square  aperture  mesh  metal  grating.  It  is  shown  that  for  a given 
QWIP  the  total  power  and  the  diffracted  angle  of  the  higher  order  TM  diffracted  waves 
depend  on  two  normalized  parameters  s = \fg  and  h = a/g,  which  are  functions  of 
wavelength  and  grating  parameters.  By  scaling  the  universal  plots  shown  in  Fig.  3.3 
and  Fig.  3.5,  the  optimum  grating  period  and  aperture  width  can  be  obtained  for 
any  specific  infrared  spectrum.  Moreover,  the  absorption  constant  and  the  coupling 
quantum  efficiency  for  a 2-D  metal  grating  coupled  QWIP  can  be  calculated  from 
these  universal  plots. 
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Figure  3.1  Schematic  diagram  for  a 2-D  square  aperture  mesh  metal  grat- 
ing coupler,  (a)  Top  view,  (b)  Side  view  under  back  illumi- 
nation. (c)  Side  view  under  front  illumination. 
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Figure  3.2  The  energy  band  diagram  for  a bound-to-miniband  transition 
GaAs/AlGaAs  QWIP.  The  ground  state  Eg  is  calculated  to 
be  290meV  after  taking  the  exchange  energy  [61,62]  into  con- 
sideration. 
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normalized  wavelength  s = ^/g 


Figure  3.3  Total  power  of  the  first  order  diffracted  waves  versus  normal- 
ized wavelength  s for  different  values  of  h.  Normal  incidence 
and  back  illumination  are  assumed. 
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normalized  wavelengths  = ^jg 


Figure  3.4  Total  power  of  the  second  order  diffracted  waves  as  a function 
of  s for  different  values  of  h.  Normal  incidence  and  back 
illumination  are  assumed. 
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normalized  wavelength  S = ^/g 


Figure  3.5  Cosine  of  the  angle  between  the  electric  field  vector  of 
diffracted  waves  in  the  quantum  well  region  and  the  normal 
of  the  quantum  well  layers  z versus  normalized  wavelength. 
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wavelength  (X,  jjm) 


Figure  3.6  Coupling  characteristic  curve  for  the  square  aperture  mesh 
metal  grating  with  (a)  g = 4.0/^m,  a = 2.0/um.,  and  (b)  g = 
3.3//m,  a = 2.0//m. 
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wavelength  X {jjn\) 


Figure  3.7(a)  Absorption  constant  versus  wavelength  for  a 2-D  square  aper- 
ture mesh  metal  grating  coupled  BTM  QWIP  for  two  different 
grating  periods  and  aperture  sizes  and  for  a 45°  incident  illu- 
mination. 
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wavelength  X {jjrr\) 


Figure  3.7(b)  Corresponding  coupling  quantum  efficiency  for  the  BTM 
QWIP  shown  in  Fig.  3.7(a). 


CHAPTER  4 

GEOMETRY  CONSIDERATIONS  OF  TWO-DIMENSIONAL 
GRATING  COUPLER  FOR  GaAs  QUANTUM  WELL 
INFRARED  PHOTODETECTOR 

4.1  Introduction 

Various  new  types  of  photodetector  have  been  demonstrated  by  GaAs/AlGaAs 
doped  quantum  well  [63-65]  structure,  which  rely  on  photodetection  involving  elec- 
tronic transitions  between  two  different  quantum  states  of  the  quantum  wells.  The 
operation  of  these  detectors  based  on  so  called  ‘intersubband’  absorption  of  infrared 
radiation  is  preferable  to  other  types  of  detectors  such  as  narrow  band  gap  or  Schot- 
tky  barrier  detectors  in  the  wavelength  range  of  8-12  /xm.  Due  to  the  excellent  large 
area  uniformity  [66]  of  GaAs/AlGaAs,  multiple  quantum  well  infrared  photodetec- 
tors (QWIPs)  show  a significant  promise  in  the  manufacture  of  2-D  high  performance 
focal  plane  arrays  (FPAs). 

In  an  intersubband  absorption  QWIP,  the  electrons  are  excited  from  bound 
ground  state  to  a higher  energy  level  miniband  [19,32]  or  continuous  band  [67,68], 
in  which  the  electrons  can  move  under  an  applied  electric  field  perpendicular  to  the 
plane  of  the  quantum  wells  and  thus  produce  photocurrent.  The  corresponding  tran- 
sitions are  allowed  only  when  one  of  the  electric  field  components  of  incident  IR 
radiation  is  polarized  in  the  direction  of  the  electron  motion  [69-72].  The  IR  light 
normally  irradiates  on  the  surface  of  detectors  will  not  be  absorbed  for  n-type  QWIPs. 
Therefore,  some  special  coupling  techniques  have  to  be  employed,  such  as  45°  angle 
polished  facet  on  one  edge  of  the  detectors,  Brewster’s  angle  illumination,  1-D  linear 
grating  [21,24]  and  2-D  cross  grating  [25,26].  Input  coupling  through  an  angle  pol- 
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ished  facet  increases  the  absorption  via  multiple-pa^s  total  internal  reflection  [42],  but 
severely  limits  device  geometry  for  detector  applications.  For  a simple  Brewster’s  an- 
gle incidence  a relatively  large  number  of  quantum  wells  is  required  in  order  to  achieve 
efficient  absorption.  The  formation  of  FPA  is  possible  in  1-D  linear  grating  coupler, 
but  the  coupling  is  polarization  dependent.  To  improve  the  coupling  efficiency,  we 
use  a planar  square  aperture  mesh  metal  grating.  In  our  design,  the  top  contact  sur- 
face of  the  detector  is  covered  with  a conducting  screen  perforated  periodically  with 
square  apertures  formed  by  standard  photolithography  and  lift-off  techniques.  The 
grating  scatters  the  incident  IR  light  into  new  electromagnetic  waves  which  obey  the 
absorption  rules  as  described  above. 

However,  the  fabrication  of  extremely  small  grating  structure  inevitably  involves 
cell-to-cell  photoresponse  nonuniformity  of  a FPAs,  i.e.,  individual  detector  may  ex- 
hibit differences  in  spectral  response,  which  in  turn  could  introduce  spatial  noise  [72] 
and  limit  the  sensitivity  of  IR  imaging  camera.  Performance  models  for  infrared  detec- 
tors such  as  detectivity  (D*)  or  responsivity  (Ri)  are  based  on  the  implicit  assumption 
that  all  the  detectors  in  an  IR  image  system  are  exactly  the  same.  Unfortunately, 
these  models  fail  to  explain  the  dominant  one— spatial  noise  in  high  background  envi- 
ronments (8  to  12  /im  image  of  practical  IR  staring  sensors).  As  a result,  the  predict 
performance  of  system  sensitivity  is  usually  optimistic.  Since  the  size  of  the  square 
apertures  in  the  metal  grating  is  around  one  micrometers,  some  variation  in  the  width 
and  shape  of  the  grating  aperture  pattern  can  be  expected.  The  effect  on  grating  per- 
formance by  the  variations  in  aperture  width  and  shape  must  be  considered  during 
the  QWIP  FPA  design. 

In  this  chapter  we  performed  a numerical  analysis  of  the  effects  in  the  coupling 
efficiency  by  grating  size  and  shape  variation  for  a square  aperture  mesh  metal  grat- 
ing coupler  on  GaAs  QWIP.  The  coupling  stability  in  grating  aperture  dimension 
was  discussed.  Finally,  an  aperture  dimension  which  can  avoid  inducing  uncorrected 
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spatial  noise  was  also  determined. 

4.2  Theory  and  Illumination 

The  2-D  square  aperture  mesh  metal  grating  with  aperture  width  a and  period- 
icity parameter  or  grating  constant  g deposited  on  top  surface  of  a GaAs  QWIP  is 
illustrated  in  Fig.  3.1.  The  whole  grating  arrangement  looks  like  two  groups  of  equally 
spaced  parallel  metal  lines  intersect  one  another  in  rectangular  angle.  When  a plane 
wave  light  strikes  the  grating  under  normal  incidence  illumination,  the  light  is  scat- 
tered into  several  order  diffracted  waves  whose  wave  vectors  have  parallel  components 
equal  to  integer  multiples  of  the  reciprocal  lattice  vector  of  the  grating 

K = yp  (4-1) 

K = (4.2) 

where  (p,  q)  is  the  diffracted  order,  and  p,q  = — oo, . . . , —2,  — 1, 0, 1, 2, . . . , oo.  The 
total  wave  vector  of  each  order  diffracted  waves;  including  evanescent  waves,  obeys 

(Kr  + (K)"  + = (2^/^)'  (4-3) 

where  is  the  perpendicular  wave  vector  of  the  diffracted  waves  which  is  real  for 
propagating  waves  and  pure  imaginary  for  evanescent  waves,  and  A is  the  wavelength 
of  the  incident  wave.  For  back  side  illumination  as  shown  in  Fig.  3.1(b),  the  wave 
solution  in  the  quantum  well  region  consists  of  the  incident  wave  plus  a superposition 
of  reflected  diffracted  waves.  The  wave  solution  in  the  air  region  consists  of  a su- 
perposition of  transmitted  diffracted  waves.  The  incident,  reflected  and  transmitted 
waves  are  expressed,  respectively,  by 

E'  = Aooi^OOl  + Aqo2^002  (4'4) 

2 

E = ^^^^Rpqr^pqr  (4-5) 

p,q  r=l 
2 

F*  = Wr  $ 

^ J / J ^pqr'*^pqr 

p,q  r=\ 


(4.6) 
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where  are  the  orthonormal  Floquet  modes,  Aqqt,  Rpqr  and  Tpqr  are  incidence, 
reflection  and  transmission  coefficients,  and  r = 1 or  2 is  designated  as  TE  or  TM 
Floquet  modes,  respectively.  By  matching  the  tangential  electric  field  of  each  wave 
at  the  grating-air  interface  and  elaborating  the  method  of  moments  [30,59],  the  final 
equation  system  is  transformed  to  a matrix  expression  to  be  solved  numerically  with 
the  aid  of  digit  computer.  The  coefficients  of  the  diffracted  waves  are  obtained  in 
terms  of  the  power  of  incident  wave.  The  normalized  power  carried  by  (p,  q)  order 
TM  reflected  wave  is  given  by 


Vpq2  — 


\Rpq2y 


(4.7) 


cos  0pq2 

where  0pq2  is  the  angle  between  the  diffracted  ray  and  the  grating  normal.  Figure  3.3 
shows  the  normalized  power  of  the  first  order  TM  diffracted  waves  plotted  in  the 
normalized  scales  X/g  and  a/g.  The  second  order  diffracted  waves  will  not  appear 
until  \/g  is  smaller  than  2.298  (using  = 3.25  for  GaAs  at  17 K).  The  coupling 
efficiency  of  the  grating  depends  not  only  on  the  power  of  each  higher  order  TM 
propagating  mode  but  also  on  the  absorption  angle  enclosed  by  the  electric  field 
vector  of  the  mode  and  the  free  carrier  motion  vector  z,  which  is  denoted  by  7 


cos  7p,2 


X/Ur 

9 


(4.8) 


where  rir  is  the  refractive  index  of  the  medium  in  which  the  wave  propagates.  Note 
that  the  power  absorption  ratio  of  quantum  wells  varies  directly  with  the  cos^  7 and 
the  diffracted  power  [32,73,74]. 

Grating  period  g remains  constant  during  the  photolithography  procedure,  the 
aperture  width  may,  however,  be  changed.  Figure  4.1  shows  the  coupling  power 
of  the  first-order  TM  diffracted  waves  varying  with  normalized  aperture  width  afg 
within  the  most  effective  coupling  regime  2.8  < Xjg  < 3.2.  The  results  clearly  show 
that  ajg  = 0.6  not  only  maximizes  the  coupling  power  but  also  is  insensitive  to  the 
variation  of  aperture  width.  This  means  that  afg  = 0.6  is  a stable  ratio  for  the 
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aperture  width  to  grating  period.  For  example,  the  grating  g = 3.3//m,  a = 2^m  used 
as  an  optimum  grating  in  Chapter  3 is  a relative  stable  design. 

Other  than  the  width  variation,  the  shape  of  the  aperture  may  also  be  distorted 
during  the  formation  of  grating  pattern  by  using  a contact-type  mask  aligner.  Fig- 
ure 4.2  illustrates  the  shape  distortion  of  the  grating  aperture.  In  the  extreme  case, 
the  square  aperture  becomes  a circular  hole  with  a radius  r equal  to  one  half  of  the 
square  aperture  width,  a.  Under  this  situation,  the  grating  consists  of  metal  film  per- 
forated with  circular  aperture.  The  characteristic  curve  of  such  a grating  is  plotted 
in  Fig.  4.3.  Almost  invariant  coupling  curve  is  obtained  for  varying  afg  = 0.6  to 
r/g  = 0.3,  which  implies  that  the  square  aperture  mesh  metal  grating  with  a normal- 
ized aperture  width  aj g = 0.6  is  insensitive  to  the  shape  distortion  even  if  grating 
pattern  is  not  properly  processed. 

Another  approach  to  avoid  the  deformation  of  grating  aperture  is  to  use  a larger 
grating  geometry  which  certainly  has  a lower  coupling  efficiency.  For  example,  a 
square  mesh  grating  coupler  with  g = 6/xm  and  a = Zgm  instead  of  = 4/^m, 
a = 2gm  can  be  chosen  to  couple  a GaAs/AlGaAs  QWIP  for  9-ll//m  detection. 
Figure  4.4  illustrates  the  characteristic  curve  of  the  first  three  order  diffracted  waves 
for  the  g = 6/xm/a  = 3fim  grating  coupler.  This  larger  scale  grating  has  larger 
aperture  size,  and  hence  is  easier  to  fabricate  within  the  same  error  tolerance  limits. 
By  using  the  equation  given  in  Equation  (3.13),  the  optical  absorption  constant  and 
coupling  quantum  efficiency  of  a miniband  transport  QWIP  (Ref  19,  Fig.  3.2)  is 
shown  in  Figs.  4.5(a)  and  (b).  In  these  figures,  two  sets  of  grating  parameters  were 
included  to  facilitate  the  comparison.  The  quantum  efficiency  rj  is  calculated  by  [51] 

>,  = ^^,(1 (4.9) 

t 

where  a,  is  the  absorption  constant  of  Rh  diffracted  mode,  I is  the  total  length  of 
quantum  wells  (in  present  case,  I = (8SA/period)x(40  periods)  = 3520A),  and  Vi  is 
the  relative  coupling  power  carried  by  the  effective  ith  order  diffracted  waves. 


60 


4.3  Conclusion 

In  conclusion,  the  difficulty  in  the  fabrication  of  small  size  square  aperture  mesh 
metal  grating  for  8-12  //m  GaAs/AlGaAs  QWIP  makes  the  consideration  in  grating 
generated  cell-to-cell  nonuniformity  becoming  important.  A detailed  analysis  of  the 
effects  of  grating  size  variation  and  shape  distortion  on  the  coupling  efficiency  for  a 
2-D  square  aperture  mesh  metal  grating  coupler  have  been  carried  out.  It  is  shown 
that  grating  generated  spatial  noise  can  be  reduced  either  by  using  a relative  stable 
aperture  width  or  a larger  grating  dimension.  The  optimum  grating  has  parameters 
of  the  ratio  a/g  = 0.6. 
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Figure  4.1  Relative  coupling  power  of  the  first  order  TM  diffracted  waves 
versus  ajg  within  the  most  significant  coupling  region  2.8  < 
s — \fg  < 3.1. 
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Figure  4.2  Top  view  of  the  distorted  grating  pattern  for  a 2-D  square 
aperture  mesh  metal  grating  coupler. 
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normalized  wavelengths  = ^jg 


Figure  4.3  The  coupling  characteristic  curve  for  the  square  aperture 
grating  with  ajg  = 0.6  and  circular  aperture  grating  with 
rjg  — 0.3,  where  r/a  = 1/2. 
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Figure  4.4  The  relative  coupling  power  of  the  1st,  2nd  and  3rd  order  TM 
diffracted  waves  versus  wavelength  for  a 2-D  square  aperture 
metal  grating  with  a = 3/im  and  g = 6//m. 
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Figure  4.5(a)  Spectral  absorption  constant  of  the  1st,  2nd  and  3rd  order  TM 
diffracted  waves  for  a miniband  transport  QWIP  (Ref.  19, 
Fig.  3.2)  coupled  by  the  grating  with  coupling  characteristic 
curve  shown  in  Fig.  4.4. 
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Figure  4.5(b)  The  coupling  quantum  efficiency  corresponding  to  Fig.  4.5(a). 


CHAPTER  5 

EXPERIMENTAL  IMPLEMENTATIONS  OF 
TWO  DIMENSIONAL  SQUARE  APERTURE 
MESH  METAL  GRATING  COUPLER 

5.1  Introduction 

Intersubband  quantum  well  detectors  have  recently  been  the  subject  of  a consid- 
erable research  effort  [74-77];  properties  important  for  many  applications  are  a wide 
spectral  response  [8],  preferably  over  the  atmospheric  window  8-12  fim  and  tunability 
of  the  peak  wavelength.  From  material  viewpoint,  this  type  of  detector  is  simpler  than 
the  narrow  band  gap  semiconductors  and  thus  has  been  widely  investigated  for  lO^m 
detection.  However,  quantum  well  absorption  takes  place  only  if  the  electric  field 
vector  of  the  IR  radiation  has  a component  perpendicular  to  the  quantum  well  layers 
[39,78,79].  This  makes  the  detector  insensitive  to  radiation  entering  in  the  direction 
normal  to  the  detector  surface  plane.  In  contrast,  a response  in  large  area  normal 
incident  illumination  is  necessary  for  the  fabrication  of  focal  plane  arrays.  Therefore, 
varieties  of  coupling  techniques  have  been  employed,  such  as  incorporating  a 1-D  lin- 
ear grating  coupler  [15],  a 2-D  metal  cross  grating  coupler  [25],  or  polishing  a 45° 
facet  in  one  edge  of  the  detector  substrate  [42].  A 2-D  double  period  planar  mesh 
metal  grating  possesses  exclusive  advantages  over  other  techniques,  which  include  (i) 
It  effectively  scatters  the  normal  incident  IR  radiation  with  transverse  electric  field 
into  a longitudinal  electric  field  to  improve  the  quantum  well  absorption  indepen- 
dent of  light  polarization,  (ii)  It  is  easy  to  fabricate  such  a planar  structure  by  metal 
deposition  on  the  QWIP  surface  and  thus  avoids  costly  etching  process. 
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In  the  grating  design,  the  square  aperture  width  ‘a’  and  grating  period  ‘p’  are 
two  key  parameters  which  govern  the  coupling  efficiency  of  the  grating.  It  has  been 
theoretically  shown  in  Chapter  4 that  the  power  coupling  for  such  a planar  grating 
structure  is  a function  of  the  ratio  of  ‘a’  to  ‘5^’,  which  is  so  called  ‘strip  factor’. 
Figure  4.1  shows  that  for  an  effective  coupling,  the  value  of  strip  factor  {h  = a/g)  for 
a square  aperture  mesh  metal  grating  must  fall  between  0.58  and  0.62.  In  this  chapter, 
a study  of  experimental  details  in  GaAs  BTM  QWIP  optically  enhanced  by  using  the 
square  aperture  mesh  metal  grating  couplers  with  parameters  g = 3.3/xm,  a — 2fim 
and  g = 4/im,  a = 2gm  under  front  and  back  side  normal  incidence  illumination  has 
been  performed. 

5.2  Sample  Preparation 

The  detector  device  sample  was  grown  by  molecular  beam  epitaxy  (MBE)  tech- 
nique. A l-//m-thick  highly  n^-doped  (1.4  x 10^*  cm“^)  GaAs  bottom  contact  layer 
was  first  grown  on  a lattice  matched  semi-insulating  GaAs  substrate.  This  is  followed 
by  the  growth  of  20  periods  of  enlarged  Ino.07Gao.93As  quantum  wells  with  a well  width 
of  lOlA  and  a dopant  density  of  1.4  x 10^®  cm“^.  The  barrier  layers  on  each  side 
of  the  quantum  well  consist  of  5 periods  of  undoped  Alo.4Gao.6As(24A)/GaAs(55A) 
superlattice  layers  which  were  grown  alternatively  with  the  InGaAs  quantum  wells. 
Finally,  a second  highly  n'*'-doped  (1.4  x 10^*  cm“^)  GaAs  cap  layer  of  0.4//m  was 
grown  on  top  of  the  QWIP  to  complete  the  structure  and  facilitate  ohmic  contact. 
The  wafer  was  polished  in  SI  GaAs  substrate  side.  Arrays  of  square  mesas  of  area 
216x216//m^  were  formed  by  etching  device  sample  through  the  active  layers  down 
to  the  bottom  contact  layer.  Following  this,  a square  ring-shaped  metal  contact  with 
9//m  width  of  AuGe/Ni/Au  was  deposited  around  the  periphery  of  the  mesa  and 
alloyed  for  ohmic  contact  formation  as  shown  Fig.  5.1. 

The  contact  lithography  utilizing  a Suss  MJB3  G-line  mask  aligner  with  320- nm 
exposure  and  commercially  produced  photomask  was  employed  for  pattern  definition. 
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The  postexposure  thickness  1.6/im  positive  polarity  photoresist  Shipley  AZ-1450J  was 
applied  according  to  manufacturer’s  specifications.  The  undiluted  Shipley  MFI  319 
was  chosen  as  the  developer.  The  first  mask  used  to  generate  the  detector  mesas 
is  transparent  except  for  the  align  marker  and  mesa  area.  After  developing  and 
postbaking,  the  sample  with  mesa  pattern  was  immersed  in  an  H2O-H3PO4-H2O2 
solution  for  mesa  etching.  The  second  mask  used  for  forming  the  detector  metal 
contact  is  opaque  except  the  top  ring-shape  contact  and  bottom  common  ground 
contact  area.  The  500A-Auo.88Geo.i2?  150A-Ni  and  850A,  99.99%-purity  Au  are  in 
turn  evaporated  direct  onto  the  contact  pattern  formed  sample.  A 3.5  minutes,  450°C 
thermal  annealing  was  finally  applied  to  the  detector  sample. 

5.3  Grating  Fabrication 

In  this  section,  we  describe  the  photolithographic  fabrication  of  the  gold  film 
grating.  The  grating  grid  is  designed  to  have  2/im  by  2//m  aperture  dimension  on  a 
3.3//m  double  periodicity.  The  photomask  is  made  of  Chromium  patterns  on  glass 
with  fifty  four  periods  of  grating  grids  in  each  detector  cell.  The  metal  grating  were 
located  within  the  interior  of  the  ring-shaped  contact  by  pure  gold  deposition  on  the 
Si-doped  GaAs  top  contact  surface  followed  by  the  metal  lift-off  technique.  High 
resolution  Shipley  System  812  photoresist  is  adopted  for  generating  grating  pattern. 
Since  the  dots  in  the  grating  pattern  are  really  tiny,  maximum  adhesion  is  required, 
which  can  be  achieved  only  on  a clean,  dry  sample  surface.  Typically,  the  wafer  was 
coated  with  photoresist  within  30  minutes  following  a dehydration  bake  of  greater 
than  125°C.  Before  priming  the  812  photoresist,  a standard  adhedsion  promoter  hex- 
amethyldisilane  (HMDS)  was  coated  on  the  wafer  by  a Headway  Research  spinner 
in  a 3750  rpm  spin  speed.  The  photoresist  was  dispensed  on  the  wafer  and  spun 
in  a 4400  rpm  speed  to  achieve  1.2/im  norminal  coating  thickness.  Eight  minutes 
softbake  was  applied  in  a convection  oven  set  at  105°C  to  reduce  the  casting  solvent 
content  and  to  promote  adhesion.  After  this  softbaking,  photoresist  film  was  exposed 
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with  the  UV  light  source  in  the  complementary  pattern  photomask.  For  optimum 
performance,  the  Shipley  MFI  321  diluted  by  deionized  water  in  the  ratio  of  three  to 
one  was  adopted  as  the  developer.  The  99.99%  purity  gold  film  of  850-A  was  then 
evaporated  onto  the  the  well-collimated  grating  pattern  by  an  electron  beam  gun. 


The  gold  film  covered  sample  wa.s  finally  immersed  in  aceton  in  an  ultrasonic  cleaner 


a GaAs  step-BTM  (SBTM)  QWIP.  The  grating  is  characterized  by  the  parameters 
of  square  aperture  width  a = 2fim  and  grating  period  g = 3.3//m.  The  IR  radia- 
tion is  incident  normally  on  the  grating  coupled  QWIP.  The  square  apertures  are 
arranged  in  a square  symmetry  geometry.  The  combination  of  normal  incidence  and 
square  symmetry  ensures  that  the  grating  coupling  are  independent  of  polarization. 
To  model  the  unknown  electromagnetic  waves  near  the  grating,  two  sets  of  orthonor- 
mal  functions,  Floquet  mode  and  waveguide  mode  n,^  spanned  by  the  square 
aperture,  were  employed  [59]  as  linear  independent  bases  as  described  in  Chapter  3 
and  4.  The  absorption  constant  a of  a QWIP  depends  on  cos^7  [80],  where  7 is  the 
absorption  angle.  In  accordance  with  this,  only  propagating  TM  Floquet  modes  of 
the  grating  scattered  field  is  eligible  for  the  intersubband  absorption  [15].  A larger 
value  of  cos^  7 implies  a larger  absorption  constant.  Besides,  the  coupling  quantum 
efficiency  of  the  grating  coupled  QWIP  depends  not  only  on  7 but  also  on  the  power 
of  effective  propagating  waves  (namely,  higher  order  TM  Floquet  modes)  which  can 
be  expressed  by 


to  stripe  the  remaining  resist,  unwanted  metal  and  leave  the  desired  grid  geometry 
delineated.  Substracting  the  area  of  contact  ring  and  bonding  pad,  the  grating  active 
area  of  the  detector  is  26400  //m^. 


5.4  Theory 


Figure  5.1  shows  the  2-D  square  aperture  mesh  metal  grating  on  top  surface  of 


(5.1) 
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where  Vi  is  the  normalized  power  of  the  ith  order  propagating  TM  mode,  a,  is  the 
corresponding  absorption  constant  and  / is  the  total  length  of  doped  QWs.  In  the 
present  sample,  I = (lOlA)  x (20  periods)  = 2020A.  The  relative  power  and  absorp- 
tion angle  of  the  first  two  order  diffracted  waves  varying  with  wavelength  A are  plotted 
in  Fig.  3.6.  The  second  order  diffracted  waves  begin  to  propagate  while  A is  smaller 
than  9.2^m  (using  refractive  index  = 3.25  for  GaAs  at  77  K).  In  the  figure,  the  ith 
order  absorption  angle  7,  is  the  quantity  used  to  calculate  the  absorption  constant  a,-. 
By  substituting  the  diffracted  power  Vi  in  Fig.  3.6  and  the  corresponding  a,  into  Eq. 
(5.1),  the  absorption  constant  is  plotted  in  Fig.  5.2.  A coupling  quantum  efficiency  of 
29%  is  obtained  at  A = 10.4//m  under  backside  normal  incident  illumination  for  the 
GaAs  BTM  QWIP  coupled  with  g — 3.3pm,  a — 2pm  square  aperture  mesh  metal 
grating.  It  is  worthy  of  noting  that  the  coupling  quantum  efficiency  can  be  further 
improved  by  choosing  a smaller  grating  period.  According  to  the  numerical  analysis  in 
Chapter  4,  the  optimum  grating  coupler  has  the  normalized  parameters  w 1,  > 
1,  and  a/p  = 0.6,  where  Ap  is  the  peak  absorption  wavelength.  Therefore,  a grating 
with  parameters  g = 3pm  and  a — 1.8pm  should  yield  a maximum  coupling  quantum 
efficiency  for  this  GaAs  SBTM  QWIP. 

5.5  Optical  Measurement 

The  photocurrent  measurement  was  carried  out  by  a globar  source,  a ORIEL 
77250  single  grating  monochromator,  a current  preamplifier  and  a chopped  lock-in 
amplifier  combination.  This  lock-in  chopper  technique  is  employed  to  enhance  the 
signal  and  block  out  the  noise  for  narrow  band  detection,  which  greatly  improve  the 
sensitivity  of  the  measurement  system.  In  addition,  the  single  grating  monochromator 
is  used  to  pass  the  desired  wavelength  of  the  IR  source.  In  order  to  eliminate  the  back- 
ground radiation,  the  IR  light  source  is  modulated  by  a chopper  at  certain  fraquency 
and  the  detected  modulated  signal  is  sent  to  the  lock-in  amplifier  with  the  chopper’s 
phase.  The  lock-in  amplifier  functions  as  both  a phase  sensitive  detector  and  a narrow 
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bandwidth  amplifier  which  amplifies  only  the  signal  at  chopper’s  frequency  and  filters 
out  the  unmodulated  signal.  Higher  frequency  will  lower  the  incident  power  accessing 
to  the  detector. 

A bias  voltage  is  required  for  the  photoconductive  operation  QWIPs.  On  the 
other  hand,  noise  is  the  limit  of  small  signal  detection  in  the  applications  of  long 
wavelength  QWIPs.  The  noise  current  from  the  TTL  power  supply  may  be  too 
large  to  dominate  the  total  noise  current.  To  solve  this,  the  battery  will  be  a much 
better  choice  for  voltage  source  because  of  its  low  noise  level  in  comparison  with  that 
of  a TTL  power  supply.  To  supply  the  desired  bias  voltage  an  operation  amplifier 
circuit  of  OPAMP  LF356  including  a variable  resistance  is  employed.  The  QWIP  was 
connected  to  the  invert  terminal,  and  the  bias  voltage  is  fed  into  the  noninvert  input 
of  the  OP  AMP  through  a voltage  divider.  The  invert  input  will  show  up  the  same 
bias  voltage  as  the  noninvert  input  because  of  the  virtual  ground  between  both  the 
OP  AMP  inputs. 

The  sample  was  mounted  in  the  cold  head  by  specially  machined  sample  holder 
inside  a liquid  helium  closed-cycle  cryogenic  system.  The  dewar  temperature  was 
measured  with  a Si  thermocouple  beneath  the  cold  head.  From  this  photocurrent 
measurement,  together  with  the  measurements  of  the  incident  light  intensity  on  the 
sample  position  by  using  a calibrated  pyroelectric  detector,  the  responsivity  varying 
with  wavelength  is  obtained. 


5.6  Results  and  Discussion 

The  intensity  of  light  source  in  long  wavelength  becomes  smaller,  as  a result, 
the  measured  photocurrent  of  the  detector  becomes  noisier.  A slight  wavelength  shift 
of  the  experimental  data  with  respect  to  the  theoretical  calculation  is  found  in  Fig. 
5.3,  since  that  the  theory  used  to  calculate  the  absorption  constant  a deals  with  the 
zero-biased  quantum  wells.  The  internal  responsivity  of  the  QWIP  Wcis  calculated 
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by  substituting  the  coupling  quantum  efficiency  given  in  Eq.  (5.1)  into  the  following 
expression 


Ri  = (5.2) 

where  the  optical  gain  go  = 0.2  was  obtained  from  the  noise  gain  measured  at  77 K 
by  the  assumptions  that  these  two  gains  are  equal.  The  bias  voltage  Vb  = —2.5V  in 
Fig.  5.3  implies  the  top  contact  is  connected  to  the  negative  polarity  of  the  battery. 
Operation  at  a higher  bias  voltage  is  unnecessary  in  order  to  avoid  a large  dark  current. 
The  grating  coupler  with  parameter  g = 3.3/xm,  a = 2/rm  shows  a larger  responsivity 
than  that  of  45°  angle  edge  detector.  The  measured  responsivity  Rp  = 490mAfW 
at  peak  response  wavelength  Ap  = 10.4/im  is  slightly  smaller  than  that  of  calculated 
value.  This  may  be  attributed  to  the  fact  that  the  conductivity  of  the  metal  grating 
(Au)  in  QWIP  device  is  smaller  than  the  theoretical  value  of  infinity.  The  current 
loss  in  conducting  mesh  metal  grating  will  certainly  lower  the  value  of  responsivity 
Rj  measured  in  real  QWIP  device. 

In  order  to  illustrate  the  effect  of  grating  period  in  the  light  coupling,  the  respon- 
sivity of  a QWIP  coupled  with  g = 4//m,  a — 2fim  grating  is  also  shown  in  Fig.  5.3. 
For  this  set  of  grating  parameters  the  coupling  characteristic  curve  is  drawn  in  Fig. 
3.6(a),  a peak  coupling  strength  occured  at  A = 12.7/im  and  fading  out  in  the  shorter 
wavelength.  Therefore  the  responsivity  of  QWIP  coupled  with  g = 3.3/xm,  a = 2^m 
is  larger  than  that  with  g = 4/im,  a = 2/zm.  The  measured  data  shown  in  Fig.  5.3 
verify  the  effect  in  responsivity  by  grating  parameters  predicted  in  Fig  3.6. 

According  to  our  previous  calculations  in  Chapter  3,  the  responsivity  of  a GaAs 
QWIP  under  back  side  illumination  is  2.34  times  higher  than  that  under  front  side 
illumination.  This  calculations  are  based  on  the  assumptions  that  the  impedance 
Zaub  looking  from  the  detector  substrate  surface  into  the  quantum  well  region  is 
equal  to  Zq/3.25  making  the  reflectivity  at  the  substrate-air  interface  is  equal  to 
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0.28.  However,  the  adding  of  the  2-D  square  aperture  metal  grating  on  top  of  the 
SBTM  QWIP  changes  Z,ub  drastically  and  exacerbates  the  light  coupling.  From  the 
experimental  results,  an  estimated  reflectivity  at  the  substrate-air  interface  is  equal  to 
0.5,  that  is,  50%  of  the  incident  power  will  be  reflected  back  to  the  air  for  a back  side 
illuminated  SBTM  QWIP  with  a 2-D  metal  grating  on  top.  As  a result,  a designed 
coating  layer  is  suggested  here  to  improve  the  back  side  light  coupling  for  this  grating 
coupled  QWIP. 


5.7  Conclusion 

In  conclusion,  we  have  performed  a experimental  study  of  the  coupling  quantum 
efficiency  and  spectral  responsivity  in  a 2-D  square  aperture  mesh  metal  grating  cou- 
pled GaAs  SBTM  QWIP  under  normal  incident  illumination.  The  coupling  efficiency 
depends  on  the  grating  period  g and  square  aperture  width  a.  The  measured  cou- 
pling quantum  efficiency  for  a GaAs  BTM  QWIP  was  found  to  be  29%  with  grating 
parameters  g = 3.3^m  and  a = 2/xm.  A further  improvement  in  the  responsivity  can 
be  achieved  by  using  a grating  with  g = 3/im  and  a = l.Sfim.  To  solve  the  strong 
reflection  in  back  shining  of  the  grating  coupled  QWIP,  a suitable  coating  material  is 
needed,  which  is  expected  to  enhance  the  responsivity  significant. 
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incident  light 


Figure  5.1  (a)  Side  view,  and  (b)  top  view  show  the  grating  diffraction, 
and  the  coordinates  of  a 2-D  square  aperture  mesh  metal  grat- 
ing coupler  formed  on  top  of  a GaAs  SBTM  detector  mesa. 


absorption  constant  (crti ) 
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Figure  5.2  Absorption  constant  versus  wavelength  for  a 2-D  square  aper- 
ture mesh  metal  grating  coupled  SBTM  QWIP  for  two  differ- 
ent grating  periods  and  aperture  sizes  and  for  a 45°  incident 
illumination.  Where  the  absorption  constant  of  the  second 
order  diffracted  waves  is  not  shown. 
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Figure  5.3  The  responsivity  of  the  2D  square  aperture  mesh  metal  grat- 
ing coupled  SBTM  QWIPs  measured  at  77  K with  an  ap- 
plied bias  voltage  Vfc,a»  = —2.51^.  The  curves  represent  (a) 
g = 3.3/im  and  a = 2/xm  under  back  side  illumination,  (b) 
g = 4//m  and  a = 2gm  under  back  side  illumination,  (c) 
g = 4//m  and  a = 2fim  under  front  illumination,  and  (d) 
launching  at  45°  polished  facet,  respectively.  The  discrete 
points  are  the  measured  data. 


CHAPTER  6 

DESIGN  OF  PLANAR  2-D  CIRCULAR  APERTURE 
METAL  GRATING  COUPLER  FOR  GaAs 
QUANTUM  WELL  INFRARED  PHOTODETECTORS 


6.1  Introduction 

Long  wavelength  infrared  detection  by  using  intersubband  optical  absorption  of 
semiconductor  quantum  well  has  been  intensively  studied  in  the  last  several  years 
[81-84].  The  excellent  large  area  uniformity  of  GaAs/AlGaAs  makes  quantum  well 
infrared  photodetectors  (QWIPs)  ideal  for  8-12  fim  long  wavelength  infrared  focal 
plane  arrays  (FPAs)  applications.  Recently,  128  by  128  large  area  preliminary  detec- 
tor arrays  for  10  /rm  infrared  imaging  have  also  been  demonstrated  [85].  However,  for 
n-type  QWIPs  the  electric  field  vector  of  the  incident  light  must  have  a component 
perpendicular  to  the  quantum  well  layers  in  order  to  induce  intersubband  absorp- 
tion [86,87].  Most  works  to  date  have  been  performed  by  shining  light  through  a 
45°  polished  facet  on  edge  of  the  detector  substrate  [42],  nevertheless,  this  coupling 
scheme  allows  only  1-D  linear  array.  For  image  applications,  planar  arrangement  of 
detectors  and  large  area  illumination  are  required.  Coupling  through  a 1-D  linear 
grating  [18,21]  on  the  QWIP  may  serve  this  purpose,  however  suffer  from  a low  cou- 
pling efficiency  due  to  its  polarization  selectivity.  A double  period  2-D  metal  grating 
[25,26,88,89]  or  random  scattering  roughened  surface  [28]  on  top  of  the  QWIP  can 
be  employed  to  deflect  the  normal  incident  light  into  an  absorbable  angle  indepen- 
dent of  light  polarization.  In  Chaps.  2 and  3 we  have  reported  the  design  of  2-D 
square  dot  reflection  metal  grating  [88]  and  2-D  square  aperture  mesh  metal  grating 
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[89]  for  n-type  QWIP.  The  unique  feature  of  these  two  metal  grating  coupler  lies  in 
the  simplicity  of  its  planar  structure,  which  can  be  fabricated  by  simple  lithography 
and  metal  deposition  techniques.  The  possible  nonuniformity  in  the  dielectric  etch- 
ing for  grating  generating  on  the  detector  surface  is  thus  eliminated,  which  greatly 
reduces  the  processing  cost.  From  a fabrication  viewpoint,  the  rectangular  corners 
of  the  square  apertures  frequently  turn  smooth  in  grating  pattern  formation,  which 
results  in  the  grating  structure  distortion.  On  the  other  hand,  no  corner  rounding 
will  be  observed  in  circular  apertures.  In  this  chapter  we  present  a numerical  and 
experimental  study  of  planar  2-D  circular  aperture  metal  grating  coupler  on  GaAs 
QWIPs.  From  the  results  shown  in  Chap.  3 that  the  grating  coupled  detector  by 
back  illumination  performs  better  than  that  by  front  illumination.  In  the  following 
discussions  we  consider  the  electromagnetic  (EM)  waves  impinging  on  a QWIP  under 
normal  incident  back  illumination. 

Light  coupling  is  polarization  independent  for  the  circular  aperture  metal  grat- 
ing in  either  square  symmetry  (Fig.  6.1(a))  or  hexagonal  symmetry  (Fig.  6.1(b)) 
arrangement.  This  is  due  to  a high  symmetry  exists  in  the  individual  unit  cell  and 
cell  arrangement  of  both  grating  structures.  It  has  been  shown  that  only  TM  com- 
ponents of  the  grating  diffracted  waves  result  in  the  intersubband  transition  in  the 
quantum  wells  [17,60].  In  addition,  the  quantum  efficiency  of  a QWIP  is  a function 
of  the  total  power  and  absorption  angle  of  the  higher-order  TM  diffracted  waves.  As 
shown  in  Fig.  6.1,  under  normal  incident  illumination,  the  normalized  total  power 
of  each  order  diffracted  waves  depends  on  the  ‘normalized  wavelength’ s = \fg  and 
the  ‘normalized  aperture  radius’  h = a/g^  where  g is  the  grating  period,  A is  the 
free  space  wavelength,  and  a is  the  radius  of  the  circular  aperture.  The  absorption 
angle  used  in  determining  the  absorption  constant  depends  on  the  parameter  ‘s’  and 
the  order  of  diffraction.  From  our  numerical  simulations,  several  universal  plots  of 
total  power  and  absorption  angle  of  higher-order  TM  diffracted  waves  versus  s and  h 
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are  developed,  which  provide  a convenient  tools  for  designing  circular  aperture  metal 
grating  coupler  on  GaAs  QWIPs. 


6.2  Basic  Theory 

The  periodicity  of  the  grating  scattering  problems  leads  to  a description  of  the 
fields  in  terms  of  the  decomposition  of  their  2-D  Floquet  mode.  Two  integers  p and 
q are  denoted  for  the  x and  y spatial  harmonics  of  the  Floquet  mode  and  a 

third  subscript  r is  used  to  denote  TE  or  TM  modes  [90,91].  A unit  power  density 
EM  wave  is  normally  incident  on  the  infinitesimally  thin  perfect  conducting  grating 
as  shown  in  Fig.  6.1(c).  The  vector  orthonormal  mode  functions  for  the  TE  and  TM 
Floquet  modes  transverse  with  respect  to  z can  be  expressed  by 


^pgi  = G(sin<^p,x  — cos  for  TE  modes 

^pg2  = G{cos(j)pqX  + sin  for  TM  modes  (6.1) 

where  the  normalized  factor  G can  be  written  in  terms  of  grating  period  g 


G 


1 

9 


for  square  symmetry 
for  hexagonal  symmetry 


(6.2) 


Besides,  all  field  quantities  are  taken  to  have  a time  varying  part  of  the  form  e-’"*, 
and  the  wave  vectors  and  in  the  x and  y axes  are  given,  respectively,  by 


k^  = A:  sin  6pg  cos  ^p,  = 
= k sin  9pg  sin  (j)pg  = 


27T 


27T 


(2?  - p) 


for  square  symmetry 

for  hexagonal  symmetry  (6.3) 
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where  6pq  and  <f>pg  are  the  corresponding  spherical  coordinates  for  the  wave  propa- 
gating vector  k of  (p,  q)-order  Floquet  mode,  and  p,q  = 0,  ±1,  ±2, . . . , ±oo.  The 
^-direction  wave  vector  depends  on  p and  q,  which  is  real  for  propagating  modes, 
pure  imaginary  for  evanescent  modes.  The  boundary  condition  needs  the  tangential 
electric  field  vectors  to  be  continuous  within  the  circular  aperture  located  at  z = 0, 
that  is 

2 2 2 

^ ^ -^OOr^OOr  “I"  ^ ^ ^ ^ Rpqr^pqr  — EE  Tpqr^pqr  withiu  the  aperture  (6.4) 

r=l  p,q  r=l  p,q  r=l 

where  the  first  summation  on  the  left-hand  side  is  the  expansion  of  the  incident  wave, 
Aoor  is  the  magnitude  of  incident  field  component  which  depends  on  the  polarization 
direction,  Rpqr  and  Tpqr  is  the  reflection  coefficient  and  transmission  coefficient  of  the 
Floquet  mode,  respectively. 

The  orthonormal  waveguide  modes  that  satisfy  the  aperture  boundary  con- 
dition  on  the  circular  aperture  itself  other  than  the  Floquet  modes  are  selected 
as  the  second  set  of  linear  independent  basis  to  represent  the  unknown  electric  field 
distribution  in  the  aperture.  As  a result,  the  transverse  electric  field  vector  is  written 
as 

EE  Tpqr^pqr  — ^ ^ ^ ^ fTmTi/flmn/ 

p,q  r m,n  I 

with  [92] 


82 


where  the  prime  superscript  represents  differential  notation,  a is  the  radius  of  the 
circular  aperture,  (p,  z)  are  the  cylindrical  coordinates  system  used,  / = 1 or  2 
denotes  TE  or  TM  mode,  respectively;  is  the  Bessel  function  of  the  first  kind  with 
order  m,  x^n  is  the  nth  zero  of  Jm{x)  and  is  the  nth  zero  of  ^Qn  ec[uals  1 

for  n = 0 and  ^on  equals  0 for  n > 1.  Assuming  that  the  impedances  looking  into  the 
quantum  well  region  from  z = 0~  plane  are  equal  to  that  of  bulk  GaAs,  the  modal 
impedances  and  rjpgr  for  TE  and  TM  waves  on  both  sides  of  the  grating  can  be 
easily  derived.  Employing  method  of  moments,  a set  of  linear  algebraic  equation  with 
the  mode  coefficients  Wmni  as  unknowns  can  be  obtained  [59] 

[M'mn,]  = |/„n4  (6.7) 

where 


yMNL 

* mnl 


r=l  'PI'' 


+ 


p,q 


Vpqr 


■)C, 


M N L*  ^mnl 


pqr 


pqr 


and 

r27T  pa 

(6.8) 

Jo  Jo 

Both  and  depend  on  grating  geometry,  and  Imni  is  the  matrix  depending 

on  incident  waves.  A principal  difficulty  in  calculating  (7^"^  by  numerical  integral 
is  the  limitation  of  CPU  time.  To  solve  this,  a closed-form  analytical  expression  of 
the  vector  function  inner  product  (7™"^  is  found  in  Appendix  A.  The  higher-order 
TM  reflection  coefficients  Rpg^  are  given  by  Eq.  (3.11).  The  normalized  power  Vpq2 
associated  with  TM  reflected  wave  diffracted  into  the  given  order  (p,  q)  is  given  by 
Eq.  (2.14).  Since  Rpgr  is  a function  of  \/g  for  a given  h,  the  normalized  power  Vpg2  is, 
also  a function  of  X/g.  The  angle  between  the  electric  field  vector  of  the  higher-order 
TM  reflected  modes  Rpg2^pq2  and  the  z direction  is  defined  as  absorption  angle  'ypq2, 
and  the  cosine  of  this  angle  depends  on  the  order  of  diffraction  and  the  normalized 
wavelength  s {—  X/g),  which  is  given  by 


cos  7p,2 


— 

Ur 


for  square  symmetry 
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s 

Ur 

where  is  the  refractive  index  of  the  detector  medium,  and  Uj.  = 3.25  for  GaAs  at 
77  K. 


v/3 


y/p^  + q'^  - 


pq 


for  hexagonal  symmetry  (6.9) 


6.3  Results  and  Discussion 

In  typical  computations,  the  lowest  52  waveguide  modes  and  5000  Floquet  modes 
with  2tt  jgyjp^  -f  q^  < 25fc  were  included.  The  addition  of  more  higher-order  modes 
made  no  noticeable  change  in  the  transmission  or  reflection  coefficients,  which  demon- 
strates the  accuracy  of  the  numerical  results.  The  problem  under  consideration  here 
is  a unit  power  density  EM  wave  impinges  normally  onto  the  2-D  metal  grating,  as 
shown  in  Fig.  6.1(c).  The  effective  coupling  of  IR  radiation  into  the  QWIP  is  due  to 
nonzeroth-order  reflected  TM  Floquet  modes,  because  the  TE  Floquet  modes  have 
an  absorption  angle  7 = 90°  and  the  evanescent  modes  produce  no  photo-signal  in 
intersubband  absorption  [17,60].  In  order  to  develop  the  generalized  relationship  of 
the  2-D  circular  aperture  metal  grating  coupler  for  QWIP,  two  normalized  grating 
parameters,  s = Xfg  and  h = a/g,  were  employed  to  facilitate  illustration. 

6.3.1  Square  Symmetry  Grating  Structure 

Figure  6.2  shows  the  first  universal  plot,  which  illustrates  the  normalized  total 
power  of  the  first-order  TM  diffracted  waves  Roi2^oi2)  ^0-12^0-125  ^-102^-102  and 
Rio2^io2  as  a function  of  s for  different  values  of  h.  The  first-order  diffracted  waves 
become  propagating  when  the  wavelength  of  the  IR  radiation  in  GaAs  is  smaller 
than  the  grating  period,  i.e.,  s = Xfg  < 3.25.  Within  the  spectral  range  in  Fig. 
6.2,  only  (0,0)  order  transmitted  waves  Tooi^oot  (*  = I52)  might  be  found  in  the 
far  field,  and  all  the  other  higher-order  transmitted  waves  were  evanescent  modes 
because  the  wavelength  in  the  air  side  is  much  greater  than  g.  In  addition,  the 
square  symmetry  arrangement  makes  the  grating  indistinguishable  between  x and  y 
directions.  The  total  normalized  power  of  the  first-order  diffracted  waves  generated 
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by  X and  y components  of  the  incident  waves  remains  the  same  for  different  incident 
polarizations.  In  other  words,  the  coupling  of  the  grating  is  polarization  independent. 
The  second-order  diffracted  waves  are  excited  by  the  grating  when  s < 2.298  with  a 
total  power  about  50%  smaller  than  that  of  the  first-order  diffracted  waves.  Figure  6.3 
illustrates  the  normalized  total  power  of  the  second-order  TM  diffracted  components 
(±1, 1)  and  (±1,  —1)  as  a function  of  s for  various  values  of  h.  A comparison  of  Figs. 

6.2  and  6.3  reveals  that  the  grating  which  excites  a larger  power  of  the  first-order 
diffracted  waves  will  prohibit  that  of  the  second-order  diffracted  waves. 

During  the  grating  fabrication,  the  grating  period  g is  certainly  invariant  while 
the  aperture  radius  a may  be  varied.  Take  this  into  consideration,  a figure  relates  the 
total  power  of  the  first-order  diffracted  waves  to  h = ajg  within  the  most  effective 
coupling  regime  2.8  < s < 3.2  is  plotted  in  Fig.  6.4.  From  the  figure,  it  is  noted 
that  the  coupling  of  grating  with  aperture  ratio  ajg  = 0.36  not  only  is  insensitive  to 
the  variation  of  aperture  radius  but  also  maximizes  the  normalized  diffracted  power 
Vejji  which  implies  that  ajg  = 0.36  is  the  optimum  aperture  radius  for  the  square 
symmetry  metal  grating. 

6.3.2  Hexagonal  Symmetry  Grating  Structure 

In  the  next  set  of  figures,  namely  Figs.  6.5,  6.6,  we  have  the  same  incident 
conditions  as  in  Figs.  6.2,  6.4  but  the  grating  has  the  hexagonal  symmetry  arrange- 
ment. Figure  6.5  shows  the  universal  plot  illustrating  the  normalized  total  power  of 
the  first-order  TM  diffracted  waves  as  a function  of  s for  different  values  of  h.  The 
six  first-order  diffracted  waves  (±1,0),  (0,±1)  and  (±1,±1)  begin  to  propagate  at 
s/3.25  = -\/3/2.  The  six  second-order  TM  diffracted  waves  (q=l,  ±1),(±1,  ±2)  and 
(±2,1)  propagate  only  for  s < 1.625  with  a total  power  about  one  tenth  of  that 
of  the  first-order  diffracted  waves.  It  should  be  noted  that  the  total  power  of  each 
order  diffracted  waves  is  independent  of  incident  polarization.  The  figure  relates  the 
total  power  of  the  first-order  diffracted  waves  to  h within  the  most  effective  coupling 
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regime  2.4  < s < 2.7  is  plotted  in  Fig.  6.6.  The  optimum  normalized  radius  for  the 
hexagonal  symmetry  metal  grating  is  obviously  a/g  = 0.31. 

6.3.3  Comparison  between  Experiment  and  Theory 

By  comparing  Figs.  6.2  and  6.5  it  is  found  that  a maximum  coupling  efficiency  of 
square  symmetry  grating  is  about  1 .08  times  larger  than  that  of  hexagonal  symmetry 
grating,  on  the  other  hand,  the  latter  has  a 1.44  times  wider  available  bandwidth  (i.e., 
the  band  between  the  normalized  wavelength  at  which  the  first  and  second  diffracted 
waves  begin  to  propagate)  than  that  of  the  former. 

The  grating  coupled  GaAs  QWIP  used  in  present  analysis  is  composed  of  20- 
period  GaAs  multiquantum  well  superlattice,  having  a well  width  of  40A,  a dopant 
density  1.2  x 10^®  cm~^,  and  a undoped  480A  barrier  of  Alo.25Gao.75As.  In  this 
QWIP,  the  peak  absorption  wavelength  is  10  /xm  and  the  photocurrent  through  the 
QWIP  is  generated  by  the  2-D  electrons  move  in  the  continuous  band  in  the  direction 
perpendicular  to  the  quantum  well  layers.  The  circular  aperture  mesh  metal  grating 
coupler  used  in  our  experiment  has  the  square  symmetry  arrangement  as  discussed  in 
Section  6.3.1.  For  effective  coupling,  the  period  of  the  grating  is  selected  so  that  the 
QWIP  absorption  band  falls  in  the  region  between  s = 2.298  and  s = 3.25  by  using 
the  first-order  TM  diffracted  waves.  Figure  6.7  shows  the  coupling  characteristic  curve 
for  the  2-D  square  symmetry  grating  with  g = 3.3/xm,  a = 1/xm,  ox  ajg  = 0.6,  which 
was  obtained  by  multiplying  the  grating  period  g — 3.3/xm  to  the  abscissa  of  curve 
ajg  = 0.6  in  Fig.  6.2.  The  cut-off  wavelength  is  10.725/xm  for  this  grating  coupler 
and  9-10.6  /xm  is  the  most  effective  coupling  regime.  The  calculated  responsivity  Rj 
for  the  QWIP  is  given  by 

ft  = - «■“')]  So  Aj 

where  Vefj  is  the  effective  coupling  power  of  the  IR  radiation;  I is  the  total  length 
of  the  doped  quantum  wells,  and  the  terms  in  the  square  brackets  is  the  expression 
of  quantum  efficiency  rj.  In  present  case,  Vgff  is  the  normalized  total  power  of  the 
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first-order  TM  diffracted  waves,  and  I = (40A/period)  x (20  periods)  = 800A.  The 
optical  gain  go  = 0.603  at  bias  Vj  = —IV  (mesa  top  negative)  is  estimated  using  the 
45°  edge-illuminated  detector.  The  intersubband  absorption  constant  is  calculated 
by  [10] 


a = cos  7 


|Mp  hukT 

ISlUr  \/Eq  + hu 


1 + ^{EF-Eo)/kT  ■ 
2 _|_ 


(6.11) 


where  7 is  the  absorption  angle  given  in  Eq.  6.9  or  Fig.  6.7,  is  the  bulk  refractive 
index,  uj  is  the  angular  frequency  of  photons,  mb  is  the  effective  mass  in  the  AlGaAs 
barriers,  Eq  is  the  ground  state  energy  level  in  the  quantum  well,  M is  the  envelop 
matrix  element,  and  a 60%  conduction  band  split  is  cissumed  [93].  The  calculated 
responsivity  for  the  grating  coupled  QWIP  with  g = 3.3/xm  and  a = 1/xm  is  shown  as 
a solid  line  in  Fig.  6.8,  and  the  measured  responsivity  corresponding  to  this  QWIP  is 
shown  by  the  square  dots  in  curve  (a)  of  Fig.  6.8.  In  this  figure,  the  grating  coupled 
QWIP  with  a/g  = l//m/3.3//m  has  a quantum  efficiency  of  17.1%  at  A = 10.2^m, 
which  is  2.7  times  higher  than  that  of  the  45°  angle  edge-illuminated  detector.  Besides, 
an  excellent  agreement  between  the  theoretical  calculations  and  the  measured  values 
of  the  responsivity  for  curve  (a)  of  Fig.  6.8.  In  order  to  facilitate  comparison,  the 
measured  responsivity  for  the  grating  coupled  QWIP  with  g = 4//m  and  a = 1.1/rm 
is  also  shown  in  Fig.  6.8  and  labeled  as  curve  (b).  The  responsivity  for  curve  (b)  is 
considerably  lower  than  that  of  curve  (a)  at  10.2/xm  wavelength.  This  is  due  to  the 
fact  that  using  grating  parameters  a/g  = lAgm/Afim  the  grating  coupling  peak  will 
shift  to  12.6  /zm  wavelength,  which  falls  beyond  the  absorption  band  of  this  QWIP. 
In  other  words,  the  grating  with  2.298  < s < 3.25  covering  the  QWIP  absorption 
band  is  a more  effective  IR  detector  coupler.  In  fact,  by  choosing  an  optimum  grating 
coupler  parameters  with  g = 3.3/zm,  a = 1.2/zm  (i.e.,  a/g  = 0.36),  further  increase  in 
the  responsivity  at  10/zm  wavelength  can  be  expected. 


6.4  Conclusion 
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In  conclusion,  we  have  performed  a numerical  simulation  and  experimental  anal- 
ysis of  spectral  responsivity  for  a 2-D  circular  aperture  mesh  metal  grating  coupled 
GaAs  QWIP.  The  total  power  and  diffracted  angle  of  the  higher-order  TM  diffracted 
waves  are  shown  to  be  related  to  two  normalized  parameters,  namely,  s = Xjg  and 
h = a/g,  which  are  functions  of  wavelength  and  grating  dimensions.  The  optimum 
grating  period  and  aperture  radius  can  be  determined  for  any  specific  infrared  spec- 
trum by  scaling  the  normalized  parameters  of  the  universal  curves  developed.  The 
responsivity  for  such  a grating  coupled  QWIP  can  also  be  calculated  from  these  uni- 
versal charts. 
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(b) 


Figure  6.1  Top  view  of  circular  apertures  mesh  metal  grating  coupler 
with  (a)  square  symmetry,  and  (b)  hexagonal  symmetry  ar- 
rangement. (c)  Side  view  showing  the  geometry  of  back  illu- 
mination. Where  g is  the  grating  period  and  a is  the  radius 
of  the  circular  aperture. 
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Figure  6.2  Normalized  total  power  of  the  first-order  diffracted  waves  of  a 
square  symmetry  grating  versus  normalized  wavelength  s for 
different  values  of  h. 
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Figure  6.3  Normalized  total  power  of  the  second-order  diffracted  waves 
of  a square  symmetry  grating  as  a function  of  s for  different 
values  of  h. 
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Figure  6.4  Total  power  of  the  first  order  diffracted  waves  versus  a/g  for 
a square  symmetry  metal  grating  coupler  at  different  values 
of  normalized  wavelength  s. 
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Figure  6.5  Normalized  total  power  of  the  first-order  diffracted  waves  of 
a hexagonal  symmetry  grating  versus  normalized  wavelength 
s for  different  values  of  h. 
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Figure  6.6  Total  power  of  the  first  order  diffracted  waves  versus  ajg  for  a 
hexagonal  symmetry  metal  grating  coupler  at  different  values 
of  normalized  wavelength  s. 
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Figure  6.7  Coupling  characteristic  curve  for  the  2-D  square  symmetry 
metal  grating  with  parameters  g = 3.3/^m,  and  o = 2.0/zm. 
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Figure  6.8  Spectral  responsivity  of  the  bound-to- continuum  QWIP  under 
various  coupling  scheme  at  bias  voltage  Vb  = —IV;  (a)  square 
symmetry  grating  with  g = 3.3/xm,  a = 1/^m,  (b)  square 
symmetry  grating  with  g = 4^m,  a = 1.1/xm,  and  (c)  45° 
polish  edge. 


CHAPTER  7 

DESIGN  OF  AN  OPTIMUM  PERFORMANCE  BTM  QWIP 

7.1  Introduction 

Recent  work  on  quantum  well  infrared  photodetectors  (QWIP)  has  demonstrated 
conclusively  that  these  detectors  can  be  implemented  in  wide  band-gap  materials  as 
a result  of  the  use  of  intraband  processes  rather  than  interband  processes  [94].  The 
quantum  well  structures  involving  such  a low-energy  transitions  utilized  for  direct 
long-wavelength  photon  detection  leads  to  both  low-voltage  operation  and  possible 
superior  performance  over  comparable  photodetectors  based  on  conventional  inter- 
band concept  [95].  The  potential  advantages  for  GaAs  QWIPs  over  conventional 
HgCdTe  detectors  in  the  long-wavelength  infrared  spectral  region  (LWIR,  8-14  fim  ) 
include  more  mature  materials  growth  and  processing  technology  with  greater  unifor- 
mity, larger  substrate  area  available  and  well-developed  passivation  skills,  possibility 
of  monolithic  integration  with  GaAs  CCDs  and  signal  processing  electronics  (and  even 
Si  electronics  with  GaAs  on  Si  growth  [83]),  better  thermal  stability,  radiation  hard- 
ness, and  inherent  high  speed  with  advantages  for  ranging.  Several  of  these  attributes 
are  particularly  desirable  for  large  area  imagine  arrays.  A detailed  understanding  in 
the  design  procedure  of  GaAs/AlGaAs  quantum  well  detectors  is,  however,  necessary, 
since  the  principle  of  this  type  of  detector  is  easy  to  extend  to  other  wide  bandgap 
materials  system  such  as  InGaAs/InAlAs  and  Si/SiGe  [82,84,96].  The  former  is  capa- 
ble of  covering  the  complete  medium  wave  infrared  band  (MWIR,  3-5  ^m).  By  using 
the  latter,  more  straightforward  integration  of  signal  processing  functions  could  be 
achieved  by  silicon-based  technology. 
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In  this  chapter,  we  describe  a theoretical  and  experimental  study  of  the  pro- 
cedures for  optimizing  the  performance  of  a BTM  QWIP.  The  effects  of  the  design 
parameters  in  QWIP  performance  including  barrier  height,  barrier  width,  quantum 
well  doping  concentration,  quantum  well  period,  and  waveguide  layer  thickness  are 
discussed.  The  results  provide  intuitive  and  quantitative  approaches  to  improve  the 
performance  of  bound-to-miniband  (BTM)  GaAs  QWIPs. 

7.2  General  Consideration  of  a BTM  QWIP 

The  ability  to  control  the  IR  spectral  bandwidth  by  varying  the  quantum  well 
parameters  is  a key  feature  for  the  QWIP  devices.  The  bound-to-continuum  (BTC) 
and  BTM  transition  QWIPs  are  the  two  QWIP  structures  which  have  been  developed 
for  the  LWIR  focal  plane  arrays  applications.  By  varying  the  well  width  and  com- 
positionally  controlled  barrier  height,  the  response  wavelengths  of  the  intersubband 
transitions  in  these  QWIPs  can  be  readily  changed  over  a wide  range  of  MWIR  and 
LWIR  spectrum.  Although  the  structures  are  quite  different,  the  operation  of  BTC 
QWIP  is  very  similar  to  that  of  a BTM  QWIP  [79].  The  infrared  radiation  is  absorbed 
in  the  doped  quantum  wells,  exciting  an  electron  from  the  highly  populated  ground 
state  and  transporting  via  the  miniband  or  continuum  states  until  it  is  collected  at 
the  electrodes  or  recaptured  into  another  quantum  well.  The  main  difference  in  the 
operation  of  these  two  types  of  detectors  is  that  the  transport  occurs  in  the  high- 
mobility  continuum  state  above  the  barriers  for  the  BTC  QWIP,  while  transport  is 
via  the  global  miniband  states  inside  the  quantum  well  for  the  BTM  QWIPs. 

The  use  of  superlattice  barrier  layer  in  a BTM  QWIP  increaises  the  barrier  height 
and  forms  a miniband  inside  the  quantum  well.  This  could  lead  to  a significant 
reduction  of  the  thermionic  emission  dark  current  over  the  barrier,  and  hence  permit 
the  raising  of  operation  temperature  for  the  QWIP.  Higher  temperature  operation  is 
desirable  for  practical  applications  of  BTM  QWIP  in  LWIR  detection  and  focal  plane 
arrays.  In  a BTM  QWIP,  the  superlattice  barriers  can  be  made  thin  to  allow  the 
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formation  of  low-lying  minibands  corresponding  to  long  wavelength,  and,  to  facilitate 
the  transport  of  carriers  through  the  superlattice.  Under  small  bias  voltage,  the 
photocurrent  is  mainly  attributed  to  the  tunneling  between  wells  and  can  be  controlled 
by  varying  the  barrier  thickness  of  the  superlattice  well,  i.e.,  thin  and  high  barriers 
allow  us  broadband  detection  with  lower  dark  current.  It  is  also  worth  noting  that 
for  photoconductive  detectors,  the  absorption  spectrum  is  found  to  be  much  broader 
for  BTC  than  for  BTM  transitions  [98,99].  As  a result,  the  absorption  constant  of 
the  former  is  significant  smaller  than  that  of  the  latter. 

7.2.1  Barrier  Height  and  AUGai_i;As  Composition 

To  reduce  the  dark  current  of  a QWIP,  a higher  barrier  height  is  needed  for  the 
quantum  well.  For  GaAs/AUGai_j;As  BTM  QWIPs,  this  can  be  achieved  by  choosing 
the  aluminum  composition  x to  be  0.38,  which  would  yield  a barrier  height  of  283 
meV  for  the  GaAs  quantum  well.  The  well  width  is  determined  so  that  there  are 
two  subband  levels  in  the  quantum  well  with  energy  separation  equals  to  the  photon 
energy  corresponding  to  an  intersubband  transition  peak  wavelength  of  around  10//m 
after  taking  the  exchange  energy  into  consideration  [61,62].  The  structure  described 
here  and  illustrated  in  Fig.  7.1  is  thus  composed  of  93A  GaAs  quantum  wells  with 
two  bound-state  located  at  Ewi,  Ew2  = 35,140  meV  above  the  conduction  band  edge 
of  GaAs  quantum  well,  respectively. 

7.2.2  Superlattice  Barrier  Width 

By  using  a short  period  superlattice  of  Alo.38Gao.62As/GaAs  as  barrier  layer  for 
the  BTM  QWIP,  the  width  of  superlattice  well  is  determined  so  that  only  one  bound- 
state  exists  in  the  superlattice  well,  which  overlaps  with  the  second  bound-state  of 
the  quantum  wells  and  creates  a global  miniband  for  carrier  transport.  A well  width 

o 

of  28 A for  the  GaAs  superlattice  well  meets  this  requirement.  The  spectral  response 
bandwidth  depends  on  the  barrier  width  of  the  superlattice  wells  [32].  For  example, 
for  a 2/zm  FWHM  absorption  bandwidth,  the  barrier  width  of  the  superlattice  well 
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is  chosen  to  be  equal  to  49A,  and  a miniband  width  of  19  meV  width  is  obtained  for 
such  a BTM  QWIP. 

A dramatic  improvement  in  the  QWIP  performance  is  achieved  by  increasing  the 
total  barrier  width  from  h = 300 A to  500 A [34].  This  is  due  to  a strong  reduction 
in  thermally  assisted  tunneling,  which  leads  to  a large  reduction  in  dark  current. 
Further  increase  of  barrier  width  gives  virtually  the  same  detectivity  D*  [28].  Since  the 
reduction  in  tunneling  due  to  the  thicker  barriers  (reduce  dark  current  Id)  is  partially 
compensated  by  the  necessity  of  using  higher  bias  voltage  (increase  Id)  to  increase 
gain.  For  this  reason,  seven  periods  Alo.38Gao,62As(49A)/GaAs(28A)  superlattice 
barrier  layer  is  chosen  as  the  optimum  value,  which  yields  a total  superlattice  barrier 
width  of  511  A. 

7.2.3  Doping  Concentration  in  the  Quantum  Well 

The  detectivity  of  a QWIP  is  nearly  constant  over  a a broad  doping  concentration 
{Nd  = 10^®  to  10^®  cm"®)  range  in  the  quantum  well  [28].  The  fabrication  of  large 
area  QWIP  arrays  can  take  advantage  of  this  doping  independent,  since  it  implies  that 
D*  will  remain  highly  uniform  across  a large  2-D  array  even  the  doping  varies  across 
the  wafer  (for  MBE  growth,  it  is  about  1%).  On  the  other  hand,  filling  of  the  charge 
storage  wells  in  the  multiplexer  circuit  by  the  dark  current  must  be  avoided,  specially 
for  large  2-D  arrays.  Therefore,  it  is  advantageous  to  lower  the  doping  concentration 
in  the  quantum  well,  which  has  essentially  the  same  detectivity  but  substantially 
lowers  dark  current.  However,  the  absorption  intensity  is  directly  proportional  to  Nd 
[32],  lowering  the  doping  will  also  reduce  the  absorption  constant,  a compromised 
doping  concentration  of  Nd  = 7 x 10^^  cm~®  is  chosen  for  optimum  performance.  In 
the  calculation  of  exchange  energy  [61,62],  a ground  state  levels  shifted  downward 
from  the  conventional  noninteracting  state  by  16  meV  was  found. 


7.3  Quantum  Well  Period 
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Based  on  the  photocarrier  gain  mechanism  [100,101],  the  photocurrent  gain  g in 
a QWIP  can  be  written  as 


1 


oc 


(7.1) 


/ JVC  + (lV+l)/i  N' 
where  L is  the  hot-carrier  recapture  mean  free  path,  I is  the  length  of  the  multiquan- 
tum well  region,  is  the  width  of  a single  quantum  well,  l\,  is  the  barrier  width  of  the 
superlattice  barrier,  and  N is  the  period  of  multiquantum  wells.  It  is  worth  noting 
that  a study  by  Stelle  et  al.  [102]  on  a set  of  standard  GaAs  QWIPs  with  4,  8,  16, 
and  32  quantum  wells  also  found  g is  inversely  proportional  to  N as  predicted  by  Eq. 
(7.1).  Increasing  N in  quantum  well  structure  produces  two  opposing  effects  (i)  the 
quantum  efficiency  r/  can  be  enhanced  by  increases  N (larger  /),  since  the  absorption 
path  increase,  and  (ii)  reducing  / can  increase  the  gain  be  increased  substantially,  since 
there  is  a decrease  in  retrapping  of  photo-electrons  into  downstream  wells.  There- 
fore it  is  useful  to  study  the  effect  of  quantum  well  period  on  the  responsivity,  dark 
current,  and  detectivity  of  BTM  QWIP.  In  this  section  we  study  the  performance 
characteristics  versus  quantum  well  period  for  three  step-BTM  QWIPs  composed  of 
20,  15,  and  4 quantum  well  periods,  labeled  20PD,  15PD  and  4PD,  respectively. 


The  20PD  multiquantum  well  structures  consist  of  lOlA  n-doped  Ino.07Gao.93As 
quantum  wells  with  doping  concentration  Nd  = 1.4  X 10^®  cm”^  and  5-period  undoped 
Alo.4Gao.6As(30  A)/GaAs  (59  A)  superlattice  barriers.  The  15PD  and  4PD  QWIPs 
were  obtained  by  etching  the  20PD  QWIP  down  to  the  sixteenth  and  fifth  well  region 
followed  by  AuGe/Ni/Au  films  deposition  to  form  ohmic  contacts.  The  dark  current 
versus  voltage  characteristics  at  77  K for  these  samples  are  shown  in  Fig.  7.2.  Since 
the  dark  current  in  these  detectors  depends  on  the  applied  electric  field  [103]  (i.e.,  the 
voltage  drop  per  period),  the  bias  voltage  H needs  to  be  converted  to  unit  of  electric 
field  to  take  into  account  the  number  of  barriers  in  each  detector.  For  this  reason, 
the  bias  voltage  for  the  20PD,  15PD  and  4PD  QWIPs  are  in  the  ratio  of  21  : 16  : 5. 

In  present  studies  the  spectral  photoresponse  was  measured  with  a globar,  a 7- 
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15/im  optical  filter,  grating  monochromator,  and  a chopper  lock-in  amplifier  system. 
The  spectral  responsivity  of  the  three  detectors,  measured  at  77  K and  the  same 
bias  field  (0.1  V/period)  are  shown  in  Fig.  7.3.  It  is  interesting  to  note  that  the 
peak  spectral  responsivity  for  the  4PD,  15PD  and  20PD  QWlPs  are  almost  identical 
which  occurs  at  Ap  = 10.5//m.  The  responsivity  determined  by  comparison  with  a 
calibrated  pyroelectric  detector  were  found  to  be  i?/(4PD)  = 0.38  A/W  (at  Vj,  = 0.5 
V),  i?/(15PD)  = 0.41  A/W  (at  H = 1.6  V),  and  il/(20PD)  = 0.33  A/W  (at  Vf,  = 2.1 
V).  The  noise  current  in  these  detectors  was  measured  at  VJ,  = 0.1  V/period.  It  is 
shown  that  in  is  limited  by  the  shot  noise  of  the  dark  current  given  by  [100] 

in  = y/ielngAf  (7.2) 


where  A/  is  the  noise  bandwidth.  Combining  the  noise  measurement  and  the  dark 
current  shown  in  Fig.  7.2  allows  the  experimental  determination  of  gain  as  ^(dPD)  = 
1.211,  ^(15PD)  = 0.196  and  5^(20PD)  = 0.103  at  VJ,  = 0.1  V/period,  respectively.  To 
relate  the  responsivity  Rj  to  quantum  well  period  N,  Rj  is  written  by 

A 1 - 

^ ~ ^^NL  + {N+l)h 


where  A:i  is  a constant  of  N,  a is  the  absorption  constant,  the  numerator  is  pro- 
portional to  quantum  efficiency  ?/,  and  the  denominator  is  equal  to  the  length  of 
the  multiquantum  well  region.  The  absorption  constant  is  determined  by  measured 
responsivity,  gain  and  those  equations  given  in  Reference  32.  The  calculated  respon- 
sivity versus  quantum  well  period  N curves  for  the  three  step-BTM  QWIPs  is  shown 
in  Fig.  7.4.  In  this  figure,  the  experimental  data  are  labeled  by  discrete  points; 
dashed  line  are  the  calculated  curves.  Similar  to  Eq.  (7.3),  the  detectivity  D*  can  be 
expressed  as 


1 — e 


al^N 


(7.4) 


>.  ■ \/Nl^  + (N  + \)k 

where  k2  is  constant  of  N and  Ad  is  the  detector  area.  The  calculated  and  measured 
detectivity  versus  periods  N are  also  shown  in  Fig.  7.4.  The  differences  between 
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measured  data  and  theoretical  prediction  can  attributed  to  the  uncertainty  in  etching 
process  in  the  4PD  and  15PD  formation.  The  exact  etching  depth  is  very  difficultly 
to  be  determined  by  using  the  simple  DEKTAK  flA  depth  profile  measurements, 
which  has  an  error  range  about  ±300A  corresponding  to  one  quantum  well  period. 
Nevertheless,  the  theoretical  predictions  has  the  same  trend  as  those  obtained  directly 
from  the  measurements.  Although  the  responsivity  of  the  4PD  QWIP  is  larger  than 
that  of  the  20PD  QWIP,  the  noise  in  Eq.  (7.2)  is  increased  by  the  square  root  of  the 
gain  y/g  and  the  detectivity  D*  is  reduced  by  this  same  factor.  As  shown  in  Fig.  7.4, 
the  optimization  periods  of  this  step-BTM  QWIP  is  20  periods,  further  increasing 
the  period  N has  little  effect  on  D*.  By  this  same  procedure,  the  optimum  period  of 
the  QWIP  design  in  Section  7.2  is  found  to  be  20. 


7.4  Blocking  Barrier  and  Dark  Current 


Using  a 450A  Alo.i6Gao.84As  layer  adjacent  to  the  end  of  superlattice  barrier  can 
effectively  block  undesirable  tunneling  carriers  [104,105].  The  quantum  wells  and 
the  blocking  layers  are  sandwiched  between  heavily  doped  GaAs  contact  layers.  The 
doping  concentration  of  the  contact  layer  Nd  = 1.0  x 10^®  cm~^  is  a suitable  density 
for  ohmic  contact  formation.  The  blocking  layer  was  designed  so  that  the  top  of  the 
barrier  is  lower  than  the  bottom  of  the  first  excited-state,  but  much  higher  than  the 
top  of  quantum  well  ground-state  to  block  the  flow  of  electrons  (dark  current)  from 
the  ground-state  to  the  collector  contact  layer.  This  idea  of  employing  such  a blocking 
layer  in  a quantum  well  IR  detector  first  appeared  in  the  work  by  Coon  et  al.  [106]. 
The  performance  improvement  by  the  blocking  layer  can  be  seen  from  the  reduction 
of  electron  transmission  coefficient  near  the  ground-state  as  shown  in  Fig.  7.5.  The 
completed  quantum  well  active  region  of  the  optimum  designed  GaAs  BTM  QWIP 
is  thus  shown  in  Fig.  7.1.  The  dark  current  density  can  be  calculated  by  [103] 


J = 


^TTqm*kT 


r \Tr\ 

Jo 


'In 


1 + e 


(EF-E)/kT 


1 _|.  ^(EF-E-qVt)/kT 


dE 


(7.5) 
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where  Vf,  is  the  bias  voltage  across  one  period  of  the  superlattice  barrier,  E is  the 
electron  energy,  Ep  is  the  Fermi  level  and  Tr  is  the  transmission  coefficient.  The  dark 
current  for  the  BTM  QWIP  of  parameters  given  in  Fig.  7.5  with  blocking  barriers  is 
estimated  to  be  below  10”®A  for  bias  voltage  less  than  0.4  V at  T = 83  K. 

7.5  Waveguide  Structure 

Using  a 2-D  planar  metal  grating  can  greatly  enhances  the  intersubband  absorp- 
tion in  QWIP  under  normal  incident  illumination.  A further  increase  in  the  coupling 
quantum  efficiency  can  be  obtained  by  using  a slab  waveguide  geometry  QWIP.  This 
structure  is  particularly  effective  for  application  where  the  total  number  of  quantum 
wells  is  small  [101].  In  such  a structure,  the  whole  QWIP  behaves  as  an  optical 
waveguide,  as  shown  in  Fig.  7.6.  The  waves  diffracted  by  the  grating  are  confined  in 
the  transverse  direction  and  guided  along  the  longitudinal  direction  to  let  radiation 
pass  the  quantum  well  stack  several  times.  The  waveguide  is  formed  by  adding  an 
internal  reflecting  layer  beneath  the  QWIP  active  region  to  reflect  the  unabsorbed 
radiation  back  into  the  quantum  well  stack.  The  cladding  layer  and  the  top  cap 
layer  are  positioned  on  each  side  of  the  quantum  well  stack.  Aluminum  arsenide  of 
refractive  index  n — 2.76  is  an  adequate  choice  as  a cladding  layer  of  GaAs  BTM 
QWIPs. 

In  this  section  we  performed  a numerical  analysis  of  GaAs/AlGaAs  grown  in- 
frared waveguide  by  using  multi-slab  waveguide  mode  theory  [107].  For  simplicity, 
we  neglect  the  radiation  loss  by  the  metal  grating  and  treat  the  waveguide  without 
grating  on  top.  As  shown  in  Fig.  7.6,  the  index  profile  of  the  waveguide  is  assumed 
to  be  invariant  in  the  direction  of  propagation,  x,  which  can  be  written  as 

e{y,z)  = eoV?{y,z)  (7.6) 

where  cq  is  the  permittivity  of  free  space  and  n{y,z)  is  the  refractive  index  profile. 
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The  electromagnetic  field  of  the  guided  mode  in  the  waveguide  is  expressed  as 

E{x,y,z)  = E“(y,z)e'‘^‘ 

H(x,y,z)  = H\y,z)e-‘l“ 

where  the  time  dependence  was  omitted  and  ^ is  the  propagating  constant  of 

the  mode.  Due  to  the  planar  geometry  of  the  waveguide  QWIP,  the  guided  fields  are 

confined  only  to  one  direction,  say  2,  and  are  independent  of  the  direction,  y,  thus 

dE^{y,z)_dH%y,z)_ 
dy  dy 

Besides,  the  planar  waveguide  supports  two  kinds  of  guided  modes:  TE  mode  has  zero 
longitudinal  electric  field  = 0),  and  TM  mode  has  zero  longitudinal  magnetic  field 
{Hx  = 0).  The  field  solution  of  TE  modes  are  readily  obtained  by  the  field  of  TM 
modes  according  to  the  concept  of  duality, 

E—^H  H ^ -E 

e — > y,  y — > t 

However,  only  TM  modes  cause  intersubband  absorption  and  are  of  major  inter- 
est, therefore,  the  following  discussion  is  focused  on  TM  modes.  First,  we  separate 
Maxwell  equation  into  transversal  component  and  longitudinal  component  as 

V«  X E°{y,z)^-  jujyH^{y,z)  = jpx  x E^{y,z) 

Vt  X H°{y,z)  - jueE^^iy^z)  =j^x  x Ht{y,z) 

There  exists  three  field  components  H^^E^  and  E°  for  the  TM  modes.  After  a com- 
parison of  components  in  three  directions,  the  above  two  equations  can  be  expressed 
in  the  following  forms 

X : 

y ■ 


, ■ OE'O  ■ L/O 


z 


(7.7) 
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Setting  H=  Hy  and  £ = utoCrE^  we  have 

dn  .. 

where  is  free  space  wavenumber.  Both  7Y  and  £ are  scalar  wave  solution  and  obey 
the  transverse  wave  equation 


H {z)  = - nXmz) 


(7.8) 


The  general  solution  of  H and  £ are 


n{z)  = ^ 

£{z)  = 


(7.9) 


where  A and  B are  constant  of  z,  and 


Y = ^ 


(7.10) 


For  one  film,  three-layer  planar  waveguide,  we  can  solve  the  wave  equation  by  setting 
the  initial  condition 

Ho  = H(Q) 

£o  = S{0) 

at  2 = 0 and  taking  H{z)  and  £{z)/er{z)  as  continuous  quantities  in  the  z direction, 
The  result  is  a matrix  expression  as  following 


where 


Ml  = 


H \ Ho 


cos  72  ^ sin  72 

■p  sin  72  cos  72 


(7.11) 


(7.12) 
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is  the  characteristic  matrix.  Extending  this  concept  to  m-layer  waveguide  with  width 
di  of  each  layer  and 

cos  7,d,  ^ sin  7,d,- 


M.  = I • • 

si  n 7,-  di  cos  7^  di 


(7.13) 


of  each  characteristic  matrix,  where  i = 1, . . . , m — 2,  the  characteristic  equation  of 
the  QWIP  waveguide  system  can  be  obtained  by 


n 

£ 


— Mm-2\d, 


m—i  \dm-2 


•Mild, 


Ho 

£0 


For  a guided  mode,  the  electromagnetic  fields  in  the  substrate  layer  are  evanescent 
wave  and  have  the  form 

n,{z)  = 

£s{z)  = -jasAe°‘‘^ 

Similar  to  the  electromagnetic  field  in  the  cap  layer,  thus 

Hciz)  = 

£’c(z)  = 

where  a,  and  Oc  are  real  number.  The  final  matrix  expression  is  given  by 


(7.14) 


B \ / mil  mi2  \ / A 

jacB  j y m2i  m22  / \ -ja^A 

and  the  dispersion  relation  for  a multilayer  slab  waveguide  used  to  solve  the  propa- 
gating constant  ^ is 


;(acmii  + a,m22)  = m2i  - a,o;cmi2 


(7.15) 


The  active  region  of  the  BTM  QWIP  under  analysis  has  20  period  quantum  wells 
with  an  intersubband  absorption  peak  at  10  ^m  wavelength  as  shown  in  Fig.  7.1. 
The  refractive  index  of  GaAs,  Alo,38Gao,62As  and  Alo.i6Gao.84As  were  taken  to  be  3.25, 
3.13  and  2.76,  respectively  [108,109].  Three  parameters  need  to  be  determined  are 
the  thickness  of  the  two  contact  layers  and  waveguide  cladding  layer.  By  adjusting 
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the  width  of  cap  and  bottom  contact  layers,  the  position  of  maximum  electric  field 
intensity  of  the  guided  mode  in  the  waveguide  can  be  shifted.  Figure  7.7  shows  the 
electric  field  distribution  of  10  /im  guided  wave  propagating  in  the  waveguide  QWIP. 
From  the  figure,  the  cap  layer  must  be  thicker  than  the  bottom  contact  layer  to 
draw  the  unbalanced  field  distribution  to  its  maximum  field  intensity  occurring  in  the 
QWIP  active  region.  The  combination  of  1.2  GaAs  cap  layer  and  0.5  //m  bottom 
contact  layer  would  yield  a much  better  performance  than  that  of  the  combination  of 
0.5  pm  cap  layer  and  1.0  pm  bottom  contact  layer.  The  optimum  widths  of  contact 
layers  are  thus  determined. 

The  thickness  of  waveguide  cladding  layer  was  chosen  to  minimize  total  layer 
thickness  without  inducing  excessive  radiative  mode  leakage.  From  the  present  anal- 
ysis, the  cladding  layer  thickness  to  make  the  field  dropping  down  to  10  % of  its 
maximum  intensity  at  the  interface  of  quantum  well/ A1  As  cladding  layer  is  a reason- 
able thickness  to  ensure  the  residual  waveguide  loss  within  QWIP  absorption  band  to 
be  very  small.  From  Fig.  7.7,  a minimum  4 pm  cladding  layer  is  chosen  to  meet  this 
requirement.  Using  the  analysis  of  Section  7.3,  the  relative  D*  versus  quantum  well 
period  of  the  GaAs  BTM  QWIP  with  waveguide  structure  of  parameters  determined 
in  Fig.  7.5  is  shown  in  Fig.  7.8.  From  the  figure,  a detectivity  maximum  is  observed 
at  = 18  for  QWIP  with  waveguide,  while,  the  maximum  detectivity  for  QWIP 
without  waveguide  requires  N larger  than  40. 

7.6  Conclusion 

We  have  performed  a theoretical  and  experiment  study  on  the  effects  of  device 
structure  and  parameters  on  the  performance  of  the  GaAs  BTM  QWIP.  A BTM 
QWIP  has  the  advantages  over  BTC  QWIP  in  that  it  can  effectively  suppress  the 
dark  current  as  well  as  increase  the  absorption  constant  significantly  by  narrowing  the 
absorption  band.  Since  the  dark  current  is  controlled  mainly  by  the  barrier  height,  by 
adding  a selective  blocking  layer  the  dark  current  can  be  further  reduced.  The  doping 
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concentration  of  the  quantum  well  has  little  effect  on  the  QWIP  detectivity.  On  the 
other  hand,  the  detectivity  can  be  improved  by  using  a large  number  of  quantum  wells 
in  the  QWIP.  From  the  present  study,  a quantum  well  period  of  20  is  the  optimum 
period  for  the  BTM  QWIP  with  doping  concentration  of  5x  10^^  cm“^.  Adding  a 
4 ^m  AlAs  cladding  layer  beneath  the  QWIP  active  region  can  greatly  increase  the 
coupling  quantum  efficiency.  Increasing  the  width  of  the  top  contact  layer  to  1.2  fim 
and  decrease  the  bottom  contact  layer  to  0.5  fim  ensures  a maximum  field  intensity 
of  IR  radiation  traveling  in  the  QWIP  active  region.  Using  the  optimized  parameters 
described  in  this  chapter,  the  D*  of  the  BTM  QWIP  can  be  increased  by  at  least  a 
factor  of  2 over  that  of  step-BTM  QWIP  without  waveguide  design  and  parameter 
optimization. 
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Figure  7.1  The  energy  band  diagram  of  an  N period  BTM  QWIP  with 
quantum  well  width  superlattice  barrier  width  h and  dark 
current  blocking  layers. 
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Figure  7.2  Dark  current  versus  bias  field  for  the  20PD,  15PD,  and  4PD 
BTM  QWIP  of  208  //m  square  mesa  measured  at  77  K. 
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Figure  7.3  Spectral  responsivity  for  the  20PD,  15PD,  and  4PD  QWIP 
measured  at  Vf,  = —0.1  V/period. 
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Figure  7.4  A comparison  of  the  calculated  and  measured  relative  respon- 
sivity  Rj  and  detectivity  D*  versus  quantum  well  period  N for 
different  period  QWIPs  without  waveguide  geometry,  where 
the  result  of  20PD  sample  is  taken  as  a perfect  match  point 
of  theoretical  calculations  and  measurements. 
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Figure  7.5  Calculated  transmission  coefficient  T*Tr  versus  the  electron 
potential  energy  E (eV)  through  the  superlattice  barrier.  The 
QWIP  consists  of  20  periods  Nu  = 7 x 10^^  cm“^  doped 
quantum  wells  of  93  A.  The  barrier  layer  on  each  side  of 
the  quantum  well  consists  of  6 period  undoped  GaAs  (2SA) / 
Alo.38Gao.62 As  (49  A)  superlattice. 
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Figure  7.6  Side  view  of  a waveguide  geometry  QWIP  shows  the  wave 
propagating  in  the  QWIP  structure. 
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thickness  z (jjm) 


Figure  7.7  The  electric  field  distribution  of  10  /im  wavelength  guided 
wave  versus  the  z axis  for  two  sets  of  different  thickness  of 
cap  and  bottom  contact  layers. 
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Figure  7.8  The  theoretical  prediction  of  detectivity  D*  versus  quantum 
well  period  N for  the  BTM  QWIPs  with  and  without  waveg- 
uide geometry  of  parameters  given  in  Fig.  7.4. 


CHAPTER  8 
SUMMARY 


The  GaAs  QWIP  becomes  a very  promising  candidate  for  long  wavelength  in- 
frared detection  due  to  the  mature  and  highly  uniform  III-V  growth  and  processing 
technologies.  In  the  detector  image  applications,  an  effective  grating  coupler  is  needed 
for  coupling  the  normal  incident  infrared  light  into  the  n-doped  quantum  wells.  In  this 
work,  we  present  three  structures  of  2-D  planar  metal  grating  couplers  for  QWIPs. 
The  exclusive  advantages  of  such  grating  couplers  are  (a)  the  coupling  is  indepen- 
dent of  light  polarization,  (b)  it  is  easy  to  fabricate  by  simple  metal  deposition.  The 
first  structure  presented  is  the  reflection  square  dot  metal  grating  with  a coupling 
efficiency  of  20%.  To  improve  the  coupling  efficiency,  a 2-D  square  aperture  metal 
grating  is  employed,  which  offers  a higher  efficiency  up  to  70%.  That  is,  70%  nor- 
mal incident  infrared  light  can  be  transferred  to  TM  mode  for  QWIP  absorption.  In 
the  third  structure,  the  2-D  circular  aperture  metal  grating  is  introduced  to  avoid 
the  difficulties  encountered  in  grating  fabrication,  which  achieves  a 70%  high  cou- 
pling efficiency,  too.  In  addition  to  the  square  symmetry  arrangement  required  in  the 
first  two  grating  structures,  the  circular  aperture  metal  grating  provides  one  more 
hexagonal  symmetry  arrangement  with  a wider  coupling  bandwidth.  In  each  of  the 
developed  structure,  the  optical  properties  of  the  metal  grating  were  studied.  The 
influences  of  grating  shape  and  dimensions  on  the  grating  coupling  efficiency  is  also 
discussed.  The  universal  plots  based  on  two  normalized  parameters  \ jg  and  a/ g offer 
an  excellent  design  tool  for  determining  the  optimum  grating  dimensions. 

A simple  approach  by  separating  grating  coupling  and  QWIP  absorption  is  de- 
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veloped  to  analyze  the  grating  coupled  QWIPs.  The  theory  for  grating  diffraction 
relies  on  the  modal  expansion  technique  and  the  method  of  moments.  A grating  of 
infinite  extension  and  perfect  conductor  is  assumed.  The  numerical  calculations  gives 
a near-exact  description  in  the  diffracted  power  and  diffracted  angle  of  the  2-D  planar 
grating.  The  selection  rules  of  intersubband  transition  proposed  that  the  absorption 
constant  is  governed  by  the  polarization  direction  of  the  light.  Together  with  the 
the  diffracted  power  density,  the  coupling  quantum  efficiency  r)  versus  wavelength 
may  be  calculated.  The  77  is  a evaluation  of  QWIP  performance  for  thermal  imaging 
applications. 

The  detector  quantum  efficiency  can  be  further  increased  by  the  use  of  waveguide 
geometry.  The  width  of  cladding  layers  as  well  as  top  and  bottom  contact  layer 
for  a QWIP  waveguide  is  calculated.  A theoretical  and  experimental  analysis  for 
quantum  well  period  is  studied,  too.  The  optimum  period  for  a structure  specified 
QWIP  provides  a economical  way  to  improve  QWIP  detectivity.  In  addition  to  the 
absorption  enhancement  by  grating,  waveguide,  and  optimum  period,  another  critical 
issue  related  to  QWIP  performance  is  the  dark  current.  The  suppression  of  dark 
current  may  be  achieved  by  a higher  superlattice  barrier  and  a selective  blocking 
layer. 

Several  works  remain  to  be  done  in  this  area.  The  experimental  results  of  grat- 
ing coupled  QWIP  shows  a strong  reflection  at  the  detector  substrate-air  interface, 
which  suggests  that  an  appropriate  coating  film  on  the  detector  back  surface  can 
greatly  improve  the  optical  coupling.  By  thinning  the  substrate  to  several  microns  to 
change  the  impedance  in  the  incident  side  is  another  approach  for  coupling  improve- 
ment. Combining  the  dark  current  suppression  and  coupling  improvement,  the  high 
sensitivity,  large  area  QWIP  FPA  can  be  expected  for  a few  years  to  come. 
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APPENDIX  A 

COMPLETE  FORM  OF  INNER  PRODUCT  OF 
FLOQUET  MODES  AND  WAVEGUIDE  MODES 
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where  G is  the  normalized  factor  as  given  in  Eq.  (2). 
(ii)  For  the  case  of  m = 0 
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